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ABSTRACT 


Future  long  dwell  high  resolution  imagery  satellites  and  space  telescopes  will 
require  very  large  flexible  primary  mirrors.  These  large  mirrors  face  many  challenges 
including  optical  surface  imperfections,  structural  vibrations,  and  jitter.  A  flexible  mirror 
can  overcome  some  of  these  challenges  by  applying  adaptive  optics  techniques  to  correct 
mirror  deformations  and  aberrations  to  produce  image  quality  data.  This  paper  examines 
and  develops  control  techniques  to  control  a  deformable  mirror  subjected  to  a 
disturbance. 

The  experimental  portion  of  the  work  uses  discrete  time  proportional  integral 
control  with  second  order  filters  to  control  disturbances  in  a  deformable  mirror  and 
correct  aberrations  in  an  adaptive  optics  system  using  laser  light.  Using  an  adaptive 
optics  testbed  containing  two  deformable  mirrors,  two  fast  steering  mirrors,  two  wave 
front  sensors,  a  position  sensor,  and  a  combination  of  lenses  the  system  corrects  a 
simulated  dynamic  disturbances  induced  in  the  deformable  mirror.  Experiments  using 
the  described  testbed  successfully  demonstrate  wavefront  control  methods,  including  a 
combined  iterative  feedback  and  gradient  control  technique.  This  control  technique 
results  in  a  three  time  improvement  in  RMS  wavefront  error  over  the  individual 
controllers  correcting  from  a  biased  mirror  position.  Second  order  discrete  time  notch 
filters  are  also  used  to  remove  induced  low  frequency  actuator  and  sensor  noise  at  0.8  Hz, 
2  Hz  and  5  Hz.  Additionally  a  2  Hz  structural  disturbance  is  simulated  on  a 
Micromachined  Membrane  Deformable  Mirror  and  removed  using  discrete  time  notch 
filters  combined  with  a  modal  iterative  closed  loop  feedback  controller,  showing  a  36 
time  improvement  in  RMS  wavefront  error  over  the  iterative  closed  loop  feedback  alone. 
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I.  INTRODUCTION 


A,  MOTIVATION 

Future  space  based  deployable  telescopes  will  be  subject  to  non- atmospheric 
disturbances.  Jitter  and  optical  misalignment  on  a  spacecraft  can  be  caused  by 
mechanical  noise  of  the  spacecraft,  and  settling  after  maneuvers.  The  introduction  of 
optical  misalignment  and  jitter  can  reduce  the  performance  of  an  optical  system  resulting 
in  pointing  error  and  contributing  to  higher  order  aberrations.  Adaptive  optics  such  as 
tip/tilt  fast  steering  mirrors  can  be  used  to  control  jitter  in  an  optical  system. 

Large  optical  surfaces  are  also  susceptible  to  local  deviation  from  perfect 
curvature  creating  higher  order  aberrations  that  require  corrections.  Future  space  based 
optics  will  be  made  from  flexible  light  weight  materials.  In  an  attempt  to  obtain  a  more 
rigid  structure,  smaller  mirrors  may  also  be  phased  together  to  create  a  larger  segmented 
mirror.  These  materials  and  large  structures  will  be  inherently  susceptible  to  surface 
errors,  vibrations,  and  noise  caused  by  both  the  environment  and  the  spacecraft.  Adaptive 
optics  concepts  and  principles  can  correct  for  aberrations  on  the  optical  surface.  Future 
space  based  large  aperture  telescopes  will  require  robust  and  responsive  control 
techniques  to  remove  dynamic  disturbances. 

In  order  to  study  these  problems,  the  Naval  Postgraduate  School  has  incorporated 
an  adaptive  optics  testbed  in  the  existing  Spacecraft  Research  and  Design  Center 
(SRDC).  This  laboratory  has  historically  studied  attitude,  pointing,  and  control  methods 
for  fine  pointing  of  optical  satellite  payloads.  The  center  has  unique  testbeds  that 
simulate  the  spacecraft  and  optical  systems  in  space  like  conditions.  This  thesis  focuses 
on  the  development  of  adaptive  optic  control  techniques  to  reduce  structural  disturbances 
in  large  aperture  optical  payloads. 
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B, 


THESIS  OBJECTIVES 


The  focus  of  this  thesis  is  to  investigate  the  adaptive  optic  control  techniques  and 
demonstrate  them  experimentally.  The  ultimate  goal  of  the  experimental  portion  is  to 
simulate  a  dynamic  disturbance  on  a  deformable  mirror  and  remove  the  disturbances 
imparted  onto  an  incoming  laser  light  source  by  removing  aberrations  in  the  wavefront 
using  deformable  mirrors  and  a  wavefront  sensor.  Knowledge  gained  from  the 
experimental  system  will  be  used  in  follow-on  research  using  a  large  light  weight 
segmented  mirror. 

C.  THESIS  OVERVIEW 

Chapter  II  provides  a  background  on  adaptive  optics  and  adaptive  optics  controls. 
An  adaptive  optics  system  is  described  using  a  discrete  time  state  space  model. 

Chapter  III  presents  the  experimental  setup  and  the  equipment  used  including 
deformable  mirrors,  fast  steering  mirrors,  and  sensors.  The  experimental  layout  is 
explained  in  detail. 

Chapter  IV  discusses  wavefront  estimation  and  the  principles  of  a  Shack- 
Hartmann  wavefront  sensor.  This  chapter  investigates  wavefront  reconstruction 
techniques  that  are  required  for  wavefront  control.  This  includes  both  indirect  wavefront 
representation  and  direct  wavefront  representation  using  modal  and  zonal  wavefront 
estimation  techniques. 

Chapter  V  discusses  and  compares  wavefront  control  techniques.  A  traditional 
iterative  closed  loop  feedback  control  technique  is  developed  using  both  direct  and 
indirect  wavefront  estimation  methods.  An  adaptive  gradient  approach  is  discussed  using 
modal  wavefront  estimation  techniques.  A  combined  iterative  feedback  and  gradient 
feedback  controller  is  also  developed.  Additionally  a  discrete  time  notch  filter  is 
developed  to  remove  known  disturbances  from  the  deformable  mirror  surface. 
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Chapter  VI  provides  experimental  results  and  analysis  of  the  wavefront  control 
techniques.  The  control  techniques  are  evaluated  from  a  biased  position  and  subject  to  a 
known  low  frequency  disturbance. 

Chapter  VII  provides  a  summary,  conclusion  and  recommendation  for  future 

work. 
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II.  BACKGROUND 


A,  ADAPTIVE  OPTICS 

Adaptive  optics  is  a  multidisciplinary  field  combing  expertise  from  physics, 
electro  optics,  controls  engineering,  mechanical  engineering,  electrical  engineering, 
materials  science  and  chemistry.  Although  many  of  the  principles  behind  adaptive  optics 
have  been  understood  for  quite  some  time  it  hasn’t  been  until  recent  times  that  adaptive 
optics  has  become  a  common  technology.  This  is  primarily  due  to  improvements  in  other 
fields  including  high  speed  computer  processing,  micro-electro-mechanical  device 
technology,  CMOS  and  CCD  cameras,  and  improved  laser  systems  (Tyson,  2000).  The 
field  continues  to  grow  as  more  challenging  applications  of  the  field  are  found.  A  recent 
example  is  the  introduction  of  large  adaptive  optics  into  space  based  telescopes  like  the 
James  Webb  Space  Telescope.  These  new  applications  introduce  new  problems  and 
challenges  that  can  only  be  solved  through  a  multidisciplinary  approach. 

The  main  purpose  of  an  adaptive  optic  system  is  to  improve  the  capability  of  an 
optical  system  by  actively  compensating  for  aberrations  in  real-time.  An  adaptive  optics 
system  can  be  simplified  to  three  subsystems.  Figure  1 .  An  active  mirror  is  the  primary 
element  where  the  surface  can  be  changed  to  match  the  phase  of  the  aberrations.  Often 
two  types  of  active  mirrors  are  employed.  A  tip/tilt  mirror  is  used  to  correct  first  order 
aberrations  while  a  deformable  mirror  is  used  to  correct  for  higher  order  aberrations.  A 
wavefront  sensor  is  the  second  element  and  is  used  to  provide  feedback  to  the  active 
mirror  in  order  to  match  the  phase  aberrations.  A  control  computer  is  the  third 
component  combining  the  wavefront  sensor  and  active  mirror  together  by  commanding 
the  actuators  of  the  deformable  mirror. 

Historically  adaptive  optics  has  been  used  in  astronomy  to  remove  wavefront 

aberrations  introduced  by  the  Earth’s  atmosphere,  in  addition  to  correcting  the  surfaces  of 

large  telescopes.  Other  traditional  adaptive  optic  applications  have  been  in  beam  control, 

particularly  in  the  use  of  high  energy  lasers  and  laser  communications.  Improving  the 

laser  wavefront  quality  offers  improvements  for  both  applications  of  directed  energy  and 
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communications  by  improving  the  efficiency  of  the  beam,  reducing  the  laser  power 
requirements  (Tyson,  2000).  Other  applications  include  optieal  relays,  both  space  based 
and  airborne.  In  both  relay  applieations  adaptive  optics  are  used  to  compensate  for  the 
distortions  eaused  by  the  atmosphere. 


▲▲ 


Figure  1  Schematic  of  Typical  Adaptive  Optics  System 

B.  ADAPTIVE  OPTICS  CONTROLS 

In  order  to  correet  for  wavefront  aberrations  regardless  of  their  eause  a  elosed 
loop  feedbaek  eontrol  law  is  required.  One  potential  problem  with  wavefront  eontrol 
using  a  large  mirror  is  control  bandwidth  separation  from  the  natural  frequency  of  the 
structure.  The  actuator  control  system  may  excite  the  structure  while  attempting  to 
control  it.  An  additional  potential  problem  is  the  low  damping  of  a  potentially  large 
space  based  mirror.  This  may  lead  to  large  resonant  peaks  in  the  frequency  response  of 
the  strueture.  These  responses  will  eontribute  to  additional  wavefront  error  as  the  surfaee 
of  the  mirror  structure  will  be  dynamieally  ehanging. 
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Adaptive  optic  control  systems  of  deformable  mirrors  require  the  use  of  multiple 
control  loops.  The  feedback  from  the  wavefront  sensor  must  be  related  to  multiple 
actuators  on  the  deformable  mirror.  What  makes  this  control  even  more  challenging  is 
the  fact  that  the  individual  actuator  control  loops  are  coupled.  Control  algorithms  can  be 
built  around  the  discrete  time  state  space  model  developed  below. 

An  adaptive  optics  system  can  be  represented  by  a  linear  state  space  model  for  the 
discrete  time  system  in  Equations  (2.1)  and  (2.2),  where  x^is  the  current  state,  w^is  the 

current  control  input,  and  is  the  output.  The  matrices  A,  B,  C,  and  D  are  the  system 
matrices  (Chen,  1993) 


(2.1) 

y,=Cx,+Du, 

(2.2) 

Applying  the  Z-transform  to  the  above  equations  and  substituting  X(z)into  T(z)  and 

setting  x(0)  =  0  results  in  Equation  (2.6). 

z[X  (z)-x(O)]  =  AX  (^z^  +  bU  (z) 

(2.3) 

X =  (^zl -  AJ  *  zx(0)  +  (z/-z4)  *5f/(z) 

(2.4) 

Y(z)  =  C(zI-Ay' zx(0)  +  fc(zI-Ay'B  +  Dlu(z) 

(2.5) 

r(z)  =  [c(z/-Ay'B  +  Dlu(z) 

(2.6) 

The  transfer  function  is  defined  in  Equation  (2.7).  The  discrete  time  system  is 
stable  if  every  bounded  input  excites  a  bounded  output  sequence.  The  transfer  function 
G(z)  is  stable  if  every  pole  of  G(z)  lies  in  the  unit  circle  of  the  z-plane  (Chen,  1993). 


G(z) 


lid 

C/(z) 


C{zl-Ay'  B  +  D 


(2.7) 
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Applying  the  discrete  time  state  space  model  in  Equations  (2.1)  and  (2.2)  to  an 
adaptive  optics  system  the  following  model,  shown  in  Equations  (2.8)  and  (2.9),  is 
developed.  In  this  model  the  state  vector,  (j) ,  is  the  wavefront  aberration,  the  matrix  B  is 
the  influence  matrix,  the  vector  c  is  the  vector  of  actuator  control  signals,  the  matrix  E  is 
a  weighting  matrix,  the  matrix  S  is  the  sensor  operator,  and  y^.  is  the  sensor  output  vector 

(Erazier  &  Tyson,  2002).  The  weighting  matrix  is  a  constant  matrix  that  weighs  the 
importance  of  the  previous  states.  In  the  adaptive  optics  system  used  in  this  thesis  the 
weighting  matrix  is  set  to  an  identity  matrix.  Therefore,  no  coupling  or  dynamics  are 
assumed  between  the  current  state  and  the  previous  state.  This  assumption  is  appropriate 
as  the  frequency  response  of  the  deformable  mirrors  used  is  very  high. 

+  (2.8) 

(2-9) 

The  above  discrete  time  state  space  model  can  also  be  used  for  a  large  mirror. 
However,  a  larger  mirror  will  have  a  lower  frequency  response  requiring  the  dynamics  to 
be  properly  modeled.  The  system  matrices  will  need  to  be  determined  experimentally  or 
by  a  finite  element  analysis.  Additional  terms  will  also  need  to  be  added  to  include  both 
the  process  noise  and  measurement  noise.  Despite  the  differences  between  the  laboratory 
mirrors  used  in  this  thesis  and  a  future  large  scale  telescope  the  control  law  development 
is  similar. 
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III.  EXPERIMENTAL  SETUP 


A,  ADAPTIVE  OPTICS  TESTBED  OVERVIEW 


Figure  2  Adaptive  Optics  Testbed  Schematic 
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Figure  3  Adaptive  Optics  Testbed 

The  adaptive  optics  testbed  is  located  in  the  Spacecraft  Research  and  Design 
Center  -  Optical  Relay  Mirror  Lab  at  the  Naval  Postgraduate  School  in  Monterey 
California.  The  components  of  the  testbed  are  mounted  on  a  Newport  Optical  Bench, 
which  can  be  floated  to  isolate  the  components  from  external  vibrations.  The  concept  is 
to  simulate  an  optical  satellite  payload  with  deformable  mirrors.  The  adaptive  optics 
testbed  uses  a  combination  of  deformable  mirrors,  tip/tilt  fast  steering  mirrors,  Shack- 
Flartmann  wavefront  sensors,  and  position  sensing  detectors  to  improve  the  quality  of  an 
imaged  object.  The  light  from  the  object  of  interest  and  a  red  reference  laser  travel 
together  through  the  optical  components  of  the  testbed.  Aberrations  can  be  input  into  the 
system  through  additional  optical  components  or  by  a  deformable  mirror.  A  Shack- 
Hartmann  wavefront  sensor  samples  the  wavefront  to  provide  feedback  to  a  control 
algorithm  which  controls  a  deformable  mirror  to  compensate  for  the  aberrations. 

The  purpose  of  the  testbed  is  to  demonstrate  advanced  control  algorithms  that 
could  be  applied  to  an  optical  payload.  The  testbed  is  set  up  with  three  different  control 
loops.  The  first  control  loop  consists  of  a  Micromachined  Membrane  Deformable  Mirror 
(MMDM)  and  Shack-lTartmann  wavefront  sensor.  This  first  loop  represents  a  primary 
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deformable  mirror  on  a  space  telescope.  The  second  control  loop  consists  of  two  Baker 
Adaptive  Optics  Fast  Steering  Mirrors  (FSM)  and  a  Position  Sensing  Detector  (PSD). 
The  second  loop  compensates  for  optical  misalignment  and  controls  tip/tilt  aberrations 
attributed  to  jitter.  The  third  control  loop  consists  of  a  Piezoelectric  Deformable  Mirror 
(PDM)  and  a  Shack-Hartmann  wavefront  sensor.  The  final  loop  is  used  to  control  higher 
order  wavefront  aberrations. 

B,  DEFORMABLE  MIRRORS 

Two  OKO  Technologies  deformable  mirrors  are  used  in  the  experimental  setup,  a 
MMDM  and  a  PDM.  The  MMDM  is  a  membrane  mirror,  with  a  5  pm  membrane 
mounted  over  a  two  dimensional  array  of  electrodes.  By  applying  a  potential  between  the 
electrodes  and  the  membrane,  the  membrane  shape  deforms.  The  PDM  is  made  from  a 
thin  solid  plate  of  glass.  The  plate  is  bonded  to  a  two  dimensional  array  of  piezoelectric 
actuators.  By  elongating  the  piezoelectric  actuators  the  mirror  deforms  (OKO 
Technologies,  2006). 


Figure  4  Simplified  MMDM  Schematic  (After,  OKO  Technologies,  2006) 


Flexible  Face  Sheet 


Figure  5 


Simplified  PDM  Schematic  (After,  OKO  Technologies,  2006) 
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The  MMDM  and  PDM  only  actuate  in  one  direction,  from  zero  voltage  due  to 
electrostatic  forces.  Therefore,  to  achieve  bi-directional  control  over  the  mirrors  they  are 
biased  to  a  half  way  point.  By  biasing  the  mirror  the  actuators  can  push  and  pull  from  the 
biased  position.  The  MMDM  and  PDM  individual  actuators  are  actuated  by  applying  an 
8-bit  control  signal  between  0  and  255,  which  is  then  applied  as  a  voltage  through  the 
mirror  electronics.  The  MMDM  surface  deflection  is  linearly  dependent  on  the  square  of 
the  applied  voltage.  The  control  signal  an  also  be  represented  as  a  value  between  +/-!. 
The  value  zero  represents  the  biased  position,  the  values  of  +/-1  represents  the  full 
positive  or  negative  actuation.  Therefore,  given  a  value,  c,  between  +/-  1  to  a  MMDM 
actuator.  Equation  (3.1)  is  used  to  compute  the  applied  control  signal.  The  PDM 
deflection  is  linearly  dependent  on  the  voltage.  Equation  (3.2).  Therefore,  the  biased 
position  for  the  MMDM  is  at  an  applied  control  signal  of  approximately  180  and  the  bias 
position  for  the  PDM  is  at  an  applied  control  signal  of  127. 

>'„™=((c  +  1).5)‘'"(255)  (3.1) 

>'™=((c  +  1).5)(255)  (3.2) 

The  MMDM  used  in  the  experimental  setup  is  15  mm  in  diameter  and  has  37 
channels,  as  shown  in  Eigure  6.  The  mirror  is  composed  of  a  silicon  chip  mounted  over  a 
holder.  The  holder  contains  the  electrode  structure  and  the  chip  contains  a  silver  nitride 
membrane,  which  is  coated  to  form  the  mirror  surface  and  grounded.  The  technical 
details  of  the  mirror  are  provided  in  Table  1  and  the  actuator  locations  are  shown  in 
Eigure  7.  The  frequency  range  of  the  mirror  is  between  0  and  500  Hz. 


Eigure  6  MMDM,  15  mm,  37  Channels 
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Parameter 

Value 

Aperture 

Circular 

Coating 

Silver  Nitride 

Aperture  Diameter 

15  mm 

Number  of  Electrodes 

37 

Control  Voltage 

0-300  Volts 

Initial  RMS  Deviation  From  Plane 

<  .45  pm 

Maximum  Deflection  at  Center 

10  pm 

Table  1  Parameters  of  MMDM  (After,  OKO  Teehnologies,  2006) 


Figure  7  MMDM  Actuator  Locations  (From,  OKO  Technologies,  2006) 


The  PDM  used  in  the  experimental  setup  is  30  mm  in  diameter  and  has  19 
actuators,  as  shown  in  Figure  8.  The  reflective  plate  is  attached  to  the  actuator  structure 
and  coated  with  a  mirror  surface.  The  technical  details  of  the  mirror  are  shown  in  Table  2 
and  the  actuator  locations  are  shown  in  Figure  9.  The  frequency  range  of  the  mirror  is 
between  0  and  1  kHz. 
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Figure  8  PDM,  30  mm,  19  Channel 


Parameter 

Value 

Aperture 

Circular 

Aperture  Diameter 

30  mm 

Number  of  Electrodes 

19 

Control  Voltage 

0-400  Volts 

Initial  RMS  Deviation  From  Plane 

<  .1  pm 

Maximum  Deflection  at  Center 

8  pm  at  400V 

Table  2  Parameters  of  PDM  (After,  OKO  Technologies,  2006) 

oil)  oil  0  12 

O  a  0  3  0  4  O  13 

OS  0  2  O  1  OS  o  14 
O  19  0  7  0  6  O  15 

O  la  o  17  O  16 

Figure  9  PDM  Actuator  Locations  (From,  OKO  Technologies,  2006) 

C.  WAVEFRONT  SENSOR 

The  Shack-FIartmann  wavefront  sensor  used  in  the  experimental  setup  is  an  OKO 
Technologies  Shack-FIartmann  wavefront  sensor.  The  sensor  includes  a  f2-inch 
Complimentary  Metal-Oxide  Semiconductor  (CMOS)  camera  model  A601f  made  by 
Easier,  Germany.  The  camera  has  656x491  pixels  and  can  operate  at  60  frames  per 
second  (fps)  and  has  a  clear  aperture  of  3.9  mm.  The  Flartmann  mask  is  made  of  fused 
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silica  and  has  a  hexagonal  geometry,  with  an  aperture  of  3.5  mm  consisting  of  127 
subaperatures  with  a  diameter  of  100  pm  (OKO  Technologies,  2006). 


Figure  10  Shaek-FIartmann  Wavefront  Sensor,  Basler  CMOS  Camera  with  Mask 

D,  FAST  STEERING  MIRRORS 

The  two  fast  steering  mirrors  (FSMs)  were  built  by  Baker  Adaptive  Opties  and 
are  shown  in  Figure  11.  The  purpose  of  the  first  FSM  is  to  introduce  disturbances  and 
vibrations  in  the  beam.  The  second  FSM  is  then  used  to  remove  the  disturbanee  by  using 
feedback  from  the  position  sensing  device.  The  FSM  has  a  one  inch  diameter  mirror  and 
uses  voice  coils  to  actuate  the  mirror.  The  voice  coils  are  placed  orthogonally  to  drive 
the  mirror  in  the  X  and  Y  directions  (tip  and  tilt).  The  FSMs  have  a  natural  frequency  at 
approximately  230  FIz  depending  on  the  direetion  of  motion  and  a  control  of  less  than 
350  Hz  depending  on  direction  of  motion.  The  FSMs  are  controlled  using  MATLAB, 
SIMULINK  and  DSPACE  by  applying  a  voltage  between  +/-5  volts. 


Figure  1 1  Baker  One  Inch  Fast  Steering  Mirror 
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E, 


POSITION  SENSING  MODULE 


Jitter  is  measured  using  a  Position  Sensing  Detector  (PSD).  This  detector  consists 
of  both  a  position  sensing  module  and  an  amplifier.  The  detector  and  amplifier  are 
manufactured  by  On-Trak  Photonics  Inc.  and  are  models  PSM2-10  and  OT-301 
respectively.  The  position  sensing  module  is  a  packaged  silicon  position  sensing 
photodiode  chip  that  produces  an  analog  output  directly  proportional  to  the  position  of  a 
light  spot  on  the  detector.  The  PSD  measures  the  reference  laser’s  position  on  the 
position  sensing  module.  The  module  is  a  10  mm  x  10  mm  duo-lateral  silicon  chip  that 
can  measure  wavelengths  between  400  and  1100  nm.  The  typical  resolution  is  250  nm 
and  the  linearity  between  analog  output  and  position  is  within  0.3%  (On-Trak  Photonics 
Inc.,  2005). 

F.  LASER 

The  reference  laser  used  is  a  red  helium  neon  laser  made  by  JDS  Uniphase,  model 
1 137P.  The  operating  wavelength  is  632.8  nm  with  a  beam  diameter  of  0.81  mm.  The 
laser  has  low  noise  with  a  maximum  RMS  of  0.2,  and  long  term  amplitude  stability  of 
2.5%  max  drift  over  8  hours.  The  beam  pointing  stability  starting  at  25°C  is  less  than 
0.10  mrads.  The  laser  provides  a  collimated  beam  of  light  that  is  used  as  a  reference  to 
determine  aberrations  of  other  optical  components  (JDSU,  2007). 

G.  OPTICAL  COMPONENTS 

Lenses  are  used  to  manage  the  reference  beam  diameter  and  ensure  that  the 
reference  laser  beam  is  collimated.  The  lenses  used  in  the  testbed  setup  include  a  20X 
microscope  objective,  and  multiple  doublets  of  different  focal  lengths.  A  microscope 
objective  is  the  first  lens  used  as  a  beam  expander  in  the  optical  path  of  the  reference 
beam.  The  microscope  objective  in  combination  with  a  doublet  lens  expands  the  beam  to 
a  one  inch  beam.  The  doublet  lenses  are  used  to  manage  diameter  of  the  beam  as  it 
travels  through  the  testbed.  A  spatial  filter  is  not  used  as  the  laser  used  produces  a  high 
quality  beam.  Beam  splitters  are  used  to  divert  a  percentage  of  the  reference  beam  in 
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order  for  the  sensors  to  provide  measurements.  Visible/infrared  two  ineh  diameter  flat 
mirrors  are  used  to  redireet  the  beam  to  different  components  of  the  testbed. 


H.  SCIENCE  CAMERA 

A  black  and  white  CCD  science  camera  is  used  to  observe  the  object  source  after 
it  has  passed  through  the  optical  components.  A  filter  is  used  to  filter  out  the  red 
reference  laser  light  source,  leaving  only  the  object  source.  The  camera  is  model  number 
IV-BWCAM3  manufactured  by  Industrial  Vision  Source.  The  camera  has  a  1/3  inch 
410,000  pixel  CCD  and  an  electronic  shutter  that  can  operate  between  60  to  100,000  Hz 
(Industrial  Vision  Source,  2004). 

I.  COMPUTER  CONTROL 

Three  desktop  computers  are  used  to  operate  the  testbed.  The  testbed  is  broken 
into  three  individual  control  loops  operated  by  three  individual  computers.  The  first 
computer  operates  the  MMDM,  which  represents  the  first  primary  optical  surface  of  the 
telescope,  and  a  wavefront  sensor.  The  second  computer  operates  the  FSMs,  PSD,  and 
the  science  camera  and  is  responsible  for  simulating  and  removing  jitter.  The  third 
computer  operates  the  PDM  and  second  wavefront  sensor  and  is  responsible  for 
providing  additional  corrections  to  the  wavefront. 

J.  DATA  ACQUISITION 

The  deformable  mirrors  are  controlled  using  MATLAB.  MATLAB  interfaces  to 
the  deformable  mirrors  through  a  MATLAB  executable  (MEX)  .dll  developed  by  Baker 
Adaptive  Optics.  A  MEX  file  is  a  dynamically  linked  subroutine  that  is  produced  from  C 
or  Eortran  source  code  (Mathworks,  2007).  In  this  case  C  source  code  is  converted  to  a 
C-MEX  file  to  provide  an  external  interface  with  the  deformable  mirrors.  The  MEX  file 
is  used  in  conjunction  with  the  OKO  Technologies  MMDM  and  PDM  drivers.  The 
individual  mirror  actuators  can  are  addressed  through  MATEAB,  and  a  control  signal 
between  0  and  255  can  be  applied  individually. 
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The  wavefront  sensors  are  also  interfaced  with  MATLAB  but  use  a  C-MEX  .dll  to 
a  memory  mapped  file.  An  executable  file  called  BAOGrabActivate.exe,  also  developed 
by  Baker  Adaptive  Optics,  is  used  to  perform  the  continuous  image  capturing  directly  to 
the  computer  RAM  via  the  memory  mapped  file.  This  allows  MATLAB  to  interface  with 
the  Basler  A601f  camera  through  the  Busier  frame  grabber  driver  using  the  1394  firewire 
port.  Using  a  cooperative,  two  part  data  acquisition  (.dll  and  .exe)  lends  itself  to  hyper- 
threaded  and  multi-core  computer  processors  as  each  component  runs  as  an  independent 
and  different  thread,  but  shares  the  memory  at  full  speed.  The  result  is  that  the  data 
acquisition  from  the  Shack-Hartmann  wavefront  sensors  uses  only  a  small  percentage  of 
the  CPU  processor,  which  is  advantageous  when  introduced  into  a  control  loop. 

K.  CALIBRATION  AND  ALIGNMENT 

The  performance  of  the  adaptive  optics  system  is  dependent  on  the  reference 
signal.  In  this  case  the  reference  signal  is  a  planar  wavefront  produced  by  the  reference 
laser.  To  ensure  that  the  wavefront  is  planar  the  beam  is  expanded  using  the  microscope 
objective  and  collimated  with  a  lens.  The  collimation  of  the  beam  is  checked  using  a 
sheer  plate. 

The  collimated  beam  is  required  to  calibrate  the  Shack-Hartmann  wavefront 
sensors.  The  wavefront  sensors  operate  based  on  the  known  positions  of  the  lenslets  on 
the  Hartmann  mask  and  their  alignment  with  the  CCD.  To  calibrate  the  wavefront  sensor 
and  remove  any  tip/tilt  bias  due  to  the  optical  components,  a  collimated  beam  was  passed 
into  the  wavefront  sensors  and  a  reference  image  was  captured.  This  reference  image  is 
used  to  measure  the  phase  difference  from  a  planar  wave. 

The  fast  steering  mirrors  are  adjusted  such  that  their  un-biased  rest  position  allows 
the  reference  beam  and  object  beam  to  pass  through  the  optical  system  without  any  tip  or 
tilt.  A  two  inch  flat  mirror  is  then  adjusted  to  ensure  the  reference  beam  is  positioned  on 
the  center  of  the  PSD. 
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IV.  WAVEFRONT  ESTIMATION 


A,  WAVEFRONT  SENSING 

Wavefront  sensing  is  required  for  a  closed  loop  adaptive  optic  feedback  control 
algorithm.  Therefore,  the  wavefront  sensor  must  have  the  spatial  resolution  and  speed  to 
be  used  in  a  real-time  feedback  system.  This  is  even  more  important  for  a  flexible 
structure  where  there  is  a  dynamic  disturbance.  There  are  two  types  of  wavefront 
sensing;  direct  and  indirect.  In  a  direct  approach,  the  wavefront  is  determined  explicitly 
while  an  indirect  approach  never  represents  the  explicit  wavefront  but  rather  transforms 
the  sensor  data  to  a  control  signal  (Tyson,  1998). 

1,  Shack-Hartmann  Wavefront  Sensor 

The  Shack-Hartmann  wavefront  sensor  output  is  proportional  to  the  wavefront 
slope.  The  wavefront  phase  is  determined  from  knowledge  of  the  wavefront  slopes.  The 
principle  of  the  Shack-Hartmann  wavefront  sensor  is  shown  through  geometry  in  Figure 
12.  The  Shack-Hartmann  wavefront  sensor  consists  of  a  lenslet  array  in  front  of  a  CMOS 
sensor.  Each  hole  on  the  lenslet  array  acts  as  an  aperture,  and  since  the  source  light 
passing  through  each  lenslet  is  converging,  the  image  produced  on  the  sensor  is  an  array 
of  spots.  The  array  of  spots  is  directly  proportional  to  the  local  wavefront  tilt  at  each 
lenslet.  The  local  wavefront  slopes,  a^j  and  p.j,  correspond  to  the  x  and  y  direction 

respectively.  The  slope  can  be  determined  by  the  Shack-Hartmann  measurements 
corresponding  the  lateral  shifts,  and  ,  of  the  local  focal  point  on  the  sensor. 

Equations  (4.1)  and  (4.2)  describe  this  relationship  where  k  is  the  wavelength  of  the 
reference  light  source,  and /is  the  focal  length  of  the  lenses  in  the  lenslet  array  (Zhu,  Sun, 
Bartsch,  Ereeman  &  Eainman,  1999). 
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Figure  12  Shack-Hartmann  Wavefront  Sensor  (After,  Southwell,  1980) 


7.n  ^ 

= - Av,.,. 


(4.1) 


(4.2) 

The  phase  ean  be  determined  by  numerieal  teehniques  if  the  number  of  slope 
measurements,  M,  is  greater  than  the  number  of  unknown  phase  points.  This  ereates  an 
over  determined  problem  and  a  solution  ean  be  found  through  a  direet  least  squares 
method.  The  error  between  the  aetual  phase  and  the  estimated  phase  determines  the 
physical  limitation  of  the  adaptive  optics  system,  and  will  affect  the  overall  performance 
of  the  closed  loop  feedback  control. 


2.  Zernike  Polynomials 

Optical  phase  can  be  represented  as  a  two  dimensional  surface  over  the  aperture. 
Deviation  from  a  reference  surface  is  considered  the  wavefront  error.  The  reference 
surface  used  in  the  experimental  work  is  a  planar  wavefront.  To  interpret  optical  test 
results  it  is  easy  to  represent  the  wavefront  as  a  polynomial  series.  The  polynomial  series 
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is  shown  in  Equation  (4.3)  where  the  Zerntke  coefficients,  and  ,  completely 

describe  the  wavefront  up  to  the  order  specified  by  the  largest  m  and  n  (Frazier  and 
Tyson,  2004). 


1”  f  r'\  "  f  r'\ 

V  n=2  V  /  n=l  m=\  J 

where. 


(4.3) 


The  series  is  in  polar  coordinates  and  the  radius,  r,  is  normalized  to  the  unit  circle, 
,  where  R  is  the  aperture  radius.  Zernike  polynomials  are  orthogonal  over  the 

interior  of  a  unit  circle,  and  therefore  appropriate  for  optical  surfaces  with  circular 
apertures.  Zernike  polynomials  can  be  transformed  to  Cartesian  coordinates  through  the 


relationship. 


Table  3  shows  the  first  24  Zernike 


polynomial  terms  using  Cartesian  coordinates. 


21 


# 

n 

m 

Polynomial 

Term 

0 

0 

0 

1 

Piston 

1 

1 

1 

X 

X-Tilt 

2 

1 

1 

T 

Y-Tilt 

3 

1 

0 

-l  +  2(x^+/) 

Focus 

4 

2 

2 

2  2 

X  -y 

Astigmatism  plus 
defocus 

5 

2 

2 

2xy 

Astigmatism  plus 
defocus 

6 

2 

1 

-2x  +  3x(^x^  + 

Coma 

7 

2 

1 

-2_y  +  3y(xV/) 

Tilt 

8 

2 

0 

l-6(.v=+/)  +  6(.v=+/f 

Third-Order 

Spherical  and  Focus 

9 

3 

3 

-  3xy ^ 

Fifth-Order 

Aberration 

10 

3 

3 

3xV-/ 

Fifth-Order 

Aberration 

11 

3 

2 

-3x^  +  3_y^  +  4x^  ^x^  +  )  “  4_y^  (^x^  k ^ ) 

Fifth-Order 

Aberration 

12 

3 

2 

-6;vy  +  8xy  (^x^ k ^ ) 

Fifth-Order 

Aberration 

13 

3 

1 

3x-12x^x^  +y^^  +  10x(^x^ 

Fifth-Order 

Aberration 

14 

3 

1 

3y-12y(x^  +y^)  +  10y(x^  +/) 

Fifth-Order 

Aberration 

15 

3 

0 

-l  +  12(x"  +/)-30(x"  +/)"  +20(x"  +/)' 

Fifth-Order 

Aberration 

16 

4 

4 

x"*  -6x^_y^  +  y'^ 

Seventh-Order 

Aberration 

17 

4 

4 

4x^y  -  4xy^ 

Seventh-Order 

Aberration 

18 

4 

3 

-4x^  +12x_y^  +5x^  ^  +_y^^-15x_y^  ^  +  y^  ^ 

Seventh-Order 

Aberration 

19 

4 

3 

-12xV  +  4/  +15xV(-^^  +/)-5/(x^  +/) 

Seventh-Order 

Aberration 

20 

4 

2 

6x^  -6y^  -20x^  (x^  +  y^'j  +  lOy^  (x^  +  _v^)  +  15x^  (x^  (x^  ■*■ 

Seventh-Order 

Aberration 

21 

4 

2 

12xy-40xy(^x^  +y^^  +  30xy(^x^ 

Seventh-Order 

Aberration 

22 

4 

1 

-4x  +  30x(^x^  +  y^^-60x(^x^  +  +35x(^x^  +  y^) 

Seventh-Order 

Aberration 

23 

4 

1 

-4y  +  30y  (x  V  _y^ )  -  60y  (x^  +  y  ^ )  V  35y  (x^  + 

Seventh-Order 

Aberration 

24 

4 

0 

1  -20(x'  +  ;;' )  +  90(x'  +  y^f  -  140(x'  +  / )'  +  70(x'  +  )" 

Seventh-Order 

Aberration 

Table  3  24  Terms  of  Zemike  Polynomials  (From,  Wyant,  2003) 
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B,  WAVEFRONT  ESTIMATION  FROM  WAVEFRONT  SLOPE 

Methods  for  determining  the  phase  of  the  wavefront  are  deseribed  as  either  zonal 
or  modal.  The  methods  are  simply  two  different  models  used  to  deseribe  the  loeal  slope 
measurements  of  a  Shaek-Hartmann  wavefront  sensor.  A  zonal  method  estimates  a  phase 
value  in  a  loeal  zone  while  the  modal  method  is  based  on  a  eoeffieient  of  an  aperture 
funetion.  In  both  oases  least-squares  estimation  is  used  for  the  phase  reoonstruotion  and 
wavefront  estimation. 

1.  Zonal  Estimation 

If  the  wavefront  is  deseribed  in  optioal  path  distanoe  over  a  small  area  or  zone 
then  the  wavefront  is  oonsidered  zonal  (Tyson,  1998).  The  zonal  estimation  method  is 
adapted  for  a  speoifio  sensor  oonfiguration,  as  the  slope  oaloulations  depend  on  the  grid 
pattern  whioh  is  then  used  to  determine  the  phase.  For  a  Shaek-Hartmann  wavefront 
sensor  a  grid  oonfiguration  is  shown  in  Figure  13.  The  dots  represent  the  lenslet  looation 
while  the  lines  represent  the  x  and  y  slopes.  Eaoh  lenslet  of  the  Shaek-Hartmann  mask 
measures  both  the  x  and  y  slope  at  the  same  point.  The  wavefront  oan  be  determined 
from  slope  measurements  through  a  least  squares  fit  of  the  slope  to  a  model  of  the  phase 
given  at  the  grid  points  or  zone.  The  phase  oan  be  modeled  by  assuming  the  phase 
differenoe  between  two  grid  points  in  the  x  and  y  direotion  is  represented  by  the 


following  two  polynomials. 

=  Cq  +  CjX  +  (4.4) 

(l)  =  CQ+c^y  +  cy  (4.5) 

The  slope  is  oaloulated  by  taking  the  derivative  of  the  previous  equations. 

=c^+2c2X  (4.6) 

V’'=c,+2c2y  (4.7) 
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Figure  13  Southwell  Geometry,  Square  Hartmann  Mask  (After,  Southwell,  1978) 


The  Shaek-Hartmann  sensor  gives  two  slope  measurements  per  grid  point, 
enabling  the  determination  of  both  oi  and  02  in  (4.6)  and  (4.7).  The  relationship  between 
slope  and  phase  is  given  below  where  the  parameter  h=D/N  where  D  is  the  diameter  of 
the  aperture  and  N  represents  the  number  of  lenslets.  Eaeh  phase  point  represents  an 
equal  sub-region  of  the  area,  ,  in  the  aperture.  The  following  equations  are  used  for  a 
square  Hartmann  mask  (Southwell,  1980). 


2  h 

2  h 


,  where 


,  where 


i=l,N-l 

j=l,N 

i=l,N 

j=l,N-l 


(4.8) 

(4.9) 


The  Hartmann  mask  used  in  the  experimental  setup  is  a  127  lenslet  hexagonal 
mask  with  a  total  aperture  of  3.5  mm,  as  shown  in  Figure  14.  Slope  and  phase  ean  be 
related  in  similar  manner  using  the  hexagonal  mask  as  the  square  mask.  However, 
speeifie  attention  must  be  given  to  the  indexing  of  the  lenslets.  Slope  averages  and  phase 
differenees  are  ealeulated  in  the  x  direetion  are  ealeulated  by  row  using  adjaeent  lenslets. 
Slope  averages  and  phase  differenees  in  the  y  direetion  are  ealeulated  using  lenslets  that 
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are  aligned  vertically.  The  lenslets  are  indexed  left  to  right  beginning  with  the  top  left 
lenslet.  Equations  (4.10)  and  (4.11)  relate  slope  and  phase  for  the  127  lenslet  hexagonal 
Hartmann  mask. 


2 


2 


^i+k) 


,  where 


f/  =  ltoll4  ^ 
i  <  6,  k=0 
7<i<13,k=l 
14<i<21,k=2 
22  <  i  <  30,  k=3 
31<i<40,k=4 
41<i<51,k=5 
52  <  i  <  63,  k=6 
64  <  i  <  74,  k=7 
75  <  i  <  84,  k=8 
85<i<93,k=9 
94<i<101,k=10 
102<i<108,k=ll 
109<i<114,k=12y 


(4.10) 


where 


f/  =  ltol02 
i  <  7,  k=16,  p=0 
8  <  i  <  15,  k=18,  p=0 
16<i<24,k=20,  p=0 
25  <  i  <  34,  k=22,  p=0 
35<i<45,k=24,  p=0 
46<i<57,k=25,p=0 
58  <  i  <  68,  k=24,  p=l 
69  <  i  <  78,  k=22,  p=3 
79  <  i  <  87,  k=20,  p=5 
88<i<95,k=18,p=7 
^96<i<102,k=16,  p=9 


(4.11) 
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Figure  14  127  Lenslet  Flexagonal  Flartmann  Mask  Used  in  the  Adaptive  Opties 

Testbed 

A  least  squares  problem  can  be  formulated  using  the  above  equations  to  compute 
phase  from  slope.  Equation  (4.12)  represents  the  least  squares  problem  where  5  is  a 
vector  of  slopes,  is  a  vector  containing  the  unknown  phase  values,  D  is  a  matrix  that 

performs  the  adjacent  slope  averaging,  and  A  is  a  matrix  that  computes  the  phase 
difference  between  two  grid  points.  The  unknown  phase  values  can  be  determined  by 
taking  the  pseudo  inverse  of  A,  represented  hy  ,  and  pre-multiplying  the  right  and  left 
hand  side  of  Equation  (4.12)  resulting  in  Equation  (4.13). 

DS  =  Al  (4.12) 

(p  =  A^'DS  (4.13) 
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2. 


Modal  Estimation 


A  wavefront  that  is  described  by  coefficients  of  the  modes  of  a  polynomial 
expansion  over  the  pupil  is  considered  modal  (Tyson,  1998).  Using  the  slope 
measurements  from  the  Shack-Hartmann  wavefront  sensor  a  set  of  coefficients,  a^. ,  can 
be  obtained  that  lit  the  following  phase  expansion  of  orthogonal  functions.  In  this  case 
(x,y)  is  a  set  of  Zemike  polynomials.  Zernike  polynomials  are  used  as  they  are  a  set 

of  orthogonal  polynomials  over  the  unit  circle.  Equation  (4.14)  can  be  written  as  a  matrix 
where  the  individual  phase  points  that  describe  the  wavefront  are  contained  in  the  vector 
the  Zernike  coefficients  are  contained  in  vector  a,  and  matrix  Z contains  a 

matrix  of  the  Zernike  terms  evaluated  at  the  phase  points  x  and  y  shown  in  Equation 
(4.15). 

M 

=  (4.14) 

k=0 

^[x,y)  =  qZ  (4.15) 

The  M  phase  expansion  coefficients,  a  ,  is  solved  using  a  least  squares  estimation, 
by  taking  the  pseudo-inverse  of  Z  at  the  measured  phase  points  from  the  Shack- 
Hartmann  wavefront  sensor  and  premultiplying  both  sides  of  Equation  (4.15)  resulting  in 
Equation  (4.16).  This  reduces  the  numerical  complexity  of  the  wavefront  estimation, 
from  the  number  of  Shack-Hartmann  lenslets  to  the  number  of  expansion  terms  used. 
a  =  zV  (4.16) 

A  slope  model  can  be  obtained  by  differentiating  Equation  (4.14).  The  resulting 
relationship  for  the  slope  in  the  x  and  y  direction  are  shown  in  Equations  (4.17)  and 

(4.18) ,  respectively.  This  allows  the  slopes,  which  are  measured  by  the  Shack-Hartmann 
wavefront  sensor  to  be  directly  related  to  the  partial  derivatives  of  the  Zernike 
polynomials.  Equation  (4.17)  and  can  be  written  in  matrix  form,  as  shown  in  Equation 

(4.19) ,  where  5  is  a  vector  of  x  and  y  slopes  with  the  dimensions  2N  x  1  and  dZ  is  a 
Matrix  of  the  partial  derivates  of  the  Zemike  terms  evaluated  at  each  lenslet,  having  2N 
rows  and  M  columns.  The  Zemike  coefficients  can  be  found  by  solving  the  least  squares 
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problem  with  the  solution  given  in  Equation  (4.20).  The  benefit  of  this  model  phase 
estimation  is  that  it  does  not  require  a  zonal  phase  estimation  as  Equation  (4.16). 
However,  by  differentiating  Equation  (4.14),  the  piston  eomponent  of  the  phase, 
eoeffieient  ,  ean  not  be  determined.  This  is  not  of  eoneern,  as  all  the  other  terms  in  the 

expansion  have  a  zero  mean,  and  so  will  the  phase  without  the  piston  term  (Southwell, 
1980). 


k=i  dx 

(4.17) 

k=i  dy 

(4.18) 

S  =  [JZ]a 

(4.19) 

a  =  [dzf  S 

(4.20) 
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V.  CONTROL  METHODS 


A, 


INFLUENCE  MATRIX 


The  MMDM  response  is  proportional  to  the  square  of  the  voltage,  as  shown  in 
Equation  (5.1),  where  d  is  the  distance  deflected  and  V  is  the  applied  control  signal.  The 
PDM  response  is  linearly  proportional  to  the  control  signal.  Equation  (5.2).  To  verify 
this  relationship  a  voltage  is  applied  to  each  electrode  ranging  from  0  to  255  while 
maintaining  a  0  control  signal  on  all  other  electrodes.  Eor  each  control  signal  applied  the 
slope  of  the  wavefront  is  measured  and  plotted  as  a  function  of  the  control  signal  for  a 
specific  actuator,  as  shown  in  Eigure  15  and  Eigure  16.  This  verifies  the  linear 
relationships  described  in  Equations  (5.1)  and  (5.2). 


^  MMDM  ^  ^MMDM 


^PDM  ^PDM 


(5.1) 

(5.2) 


Eigure  15  Shack-Hartmann  Eenslet  64  X-slope  Response  vs.  MMDM  Actuator  1 
Control  Signal  (Eeft),  Shack-Hartmann  Eenslet  64  X-slope  Response  vs. 
MMDM  Actuator  1  Control  Signal  Squared  (Right) 
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Figure  16  Shack-FIartmann  Lenslet  80  X-slope  Response  vs.  PDM  Aetuator  5 
Control  Signal 

In  Figure  15  and  Figure  16,  the  mirror  response  saturates  at  higher  eontrol  signals. 
The  usable  eontrol  signal  range  of  the  mirror  is  approximately  240  for  the  MMDM  and 
1 90  for  the  PDM  as  the  surfaee  does  not  displaee  with  the  same  linear  relationship  as  it 
does  for  lower  eontrol  signals.  These  limits  must  be  considered  when  developing  the 
control  law  and  relating  the  mirror  deformation  to  a  control  signal. 

To  control  a  deformable  mirror,  the  wavefront  sensor  data  must  be  related  to  a 
control  signal  in  order  to  provide  feedback  control  in  the  adaptive  optic  system.  This 
relationship  is  established  by  creating  an  influence  matrix.  The  influence  matrix  allows 
one  to  relate  the  control  signal  of  an  actuator  to  the  change  in  the  shape  of  the  mirror. 
The  influence  matrix  is  created  by  setting  all  the  actuator  control  signals  to  zero  or  a 
biased  control  signal  and  applying  a  maximum  control  signal  to  each  actuator  and 
recording  the  response  of  the  wavefront.  The  influence  matrix  is  also  known  as  a  poke 
matrix  as  the  matrix  is  built  by  poking  each  individual  actuator.  The  wavefront  response 
can  be  represented  indirectly  as  the  sensor  response  or  directly  as  modal  coefficients.  The 
wavefront  responses  are  represented  as  column  vectors  for  each  control  channel,  and  the 
column  vectors  make  up  the  influence  matrix,  [R]  . 
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The  size  of  the  influenee  matrix  depends  on  the  number  of  sensor  measurements, 
N,  and  the  number  of  eontrol  ehannels,  M.  The  influenee  matrix  takes  the  form  of 
Equation  (5.3),  where  the  eolumn  veetors,  5 ,  represent  the  loeal  slopes  at  eaeh  lenslet. 
Using  the  Hartmann  mask  the  number  of  slope  measurements  equals  twiee  the  number  of 
lenslets,  N,  as  slope  is  measured  in  both  the  x  and  y  direetion.  Eaeh  slope  veetor  is 
defined  by  Equation  (5.4)  where  the  slopes  in  the  x  direetion  are  preeeded  by  the  slopes 
in  the  y  direetion  The  resulting  influenee  matrix  size  for  the  Shaek-Hartmann  sensors 
used  in  the  experimental  setup  are  254  x  37  and  254  x  19  for  the  MMDM  and  PDM 
respeetively. 

^  =  U  •••  (5-3) 

^  =  [si  •••  4  4  •••  4]^  (5-4) 

Onee  the  influenee  matrix  is  developed  a  relationship  between  eontrol  signals  and 
the  sensor  response  ean  be  determined  by  Equation  (5.5),  where  c  is  a  eontrol  signal 
eolumn  veetor,  B  is  the  influenee  matrix  and  5  is  a  eolumn  veetor  of  the  sensor  data 
representing  loeal  slope  measurements  at  the  Shaek-Hartmann  lenslet.  A  desired 
wavefront  ean  be  represented  as  the  veetor  s_  and  the  required  eontrol  signal  veetor  ean 
be  eomputed  by  taking  the  pseudo-inverse  of  the  influenee  matrix  using  the  Singular 
Value  Deeomposition  (SVD)  and  pre-multiplying  both  sides  of  Equation  (5.5)  resulting 
in  Equation  (5.6).  This  method  eomputes  a  least  squares  solution  when  determining  the 
eontrol  voltage  for  a  desired  wavefront. 


=  B.c_ 

(5.5) 

(5.6) 

The  influenee  matrix  ean  be  represented  using  a  zonal  or  modal  representation. 
Instead  of  representing  the  Shaek-Hartmann  data  as  slopes  in  the  influenee  matrix  the 
data  ean  be  represented  as  zonal  phase  points  using  the  phase  estimation  method 
deseribed  in  Equation  (4.13).  The  resulting  influenee  matrix  is  deseribed  by  Equation 
(5.7).  The  influenee  matrix  ean  also  be  represented  using  a  modal  phase  estimation  as 
deseribed  in  Equation  (4.16)  and  (4.20)  (modes  from  zonal  phase  and  modes  from  slope) 
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resulting  in  the  influenee  matrix  deseribed  by  Equation  (5.8),  where  is  a  matrix  of  the 
Zernike  terms  evaluated  at  eaeh  phase  point.  In  the  experimental  setup  the  matrix  has 

the  dimensions  N  x  M,  where  M  is  the  number  of  Zernike  terms  used,  and  the  matrix 
dZ  has  the  dimensions  2N  x  M.  In  both  the  zonal  and  modal  representation  the 
relationship  between  the  eontrol  voltage  and  measured  wavefront  are  shown  in  Equations 
(5.5)  and  (5.6).  Using  a  zonal  or  modal  phase  representation  of  the  wavefront  results  in 
the  voltage  to  wavefront  relationship  shown  in  Equations  (5.9)  and  (5.10). 


b,=a^db^ 

(5.7) 

B^=ZlA^DB^ovB^=dZ^B^ 

(5.8) 

E  = 

(5.9) 

c  =  Bla 

(5.10) 

Partieular  attention  should  be  given  to  the  eondition  number  of  the  influenee 

matrix,  k  =  ,  where  cr  represents  a  singular  value.  A  poorly  eonditioned  matrix 

cr  . 

mm 

leads  to  numerieal  instabilities  when  inverting  the  influenee  matrix  to  determine  the 
eontrol  signals  as  in  Equation  (5.6).  Table  4  shows  the  eondition  numbers  of  three 
influenee  matriees.  The  first  influenee  matrix  is  eonstrueted  from  sensor  slopes,  the 
seeond  from  zonal  phase  points,  and  the  third  from  modal  Zernike  eoeffieients.  The 
eondition  numbers  were  obtained  experimentally  using  the  MMDM,  indieating  that  the 
influenee  matrix  with  modal  eoeffieients  is  the  most  numerieally  stable. 


Influenee  Matrix 

Dimensions 

K 

B 

s 

254  X  37 

6.7x10^ 

127x37 

6.7x10^ 

B. 

21  x37 

221.8 

Table  4  Experimental  Condition  Number  of  Influenee  Matrix  for  MMDM 
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B, 


ITERATIVE  FEEDBACK  CONTROL 


The  first  controller  discussed  is  used  throughout  adaptive  optics  and  involves 
using  an  iterative  closed  loop  feedback  controller.  The  controller  is  similar  to  a  closed 
loop  proportional  discrete  time  integral  controller  where  a  new  control  signal  is  updated 
based  off  the  error  multiplied  by  a  proportional  gain.  The  error  is  computed  using  the 
sensor  data  which  is  used  to  estimate  the  wavefront  and  compute  the  estimated  residual 
wavefront  aberration.  The  residual  wavefront  aberration  is  related  to  a  control  signal 
using  the  influence  matrix.  The  control  signal  representing  the  error  is  then  multiplied  by 
a  gain.  The  block  diagram  is  shown  below  and  the  plant  is  represented  by  the  influence 
matrix  B,  which  models  both  the  deformable  mirror  and  the  wavefront  sensor.  This 
control  law  is  implemented  using  different  direct  wavefront  estimation  techniques 
previously  discussed  as  well  as  using  indirect  wavefront  representation. 


Figure  17  Iterative  Feedback  Control 

1.  Indirect  Iterative  Feedback  Control 

An  indirect  control  method  in  adaptive  optics  avoids  explicitly  determining  the 
wavefront.  This  can  reduce  the  number  of  numerical  steps  in  the  control  algorithm,  but 
may  create  a  poorly  conditioned  numerical  problem.  The  slope  data  measured  by  the 
Shack-Fiartmann  wavefront  sensor  can  be  directly  implemented  into  a  feedback  control 
algorithm.  A  simple  feedback  control  algorithm  is  presented  in  Equation  (5.11).  This 
algorithm  is  an  iterative  feedback  control  loop  that  constantly  updates  the  control  signal 
vector,  c ,  based  on  a  slope  measurement  vector,  ^ ,  and  the  influence  matrix,  .  The 
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variable  g  is  a  gain  between  zero  and  one.  The  influenee  matrix  is  eonstrueted  using  a 
biased  rest  voltage  for  the  aetuators  not  being  poked.  The  influenee  matrix  is  inverted 
using  a  pseudo-inverse  SVD  approaeh. 

(5.11) 


2.  Indirect  Iterative  Feedback  Control  with  Singularity  Robust  Inverse 

The  algorithm  in  Equation  (5.11)  ean  experienee  numerieal  instabilities  when 
inverting  the  influenee  matrix  .  If  the  rank  of  (matrix  size  2N  x  M)  is  less  than  the 

minimum  of  [2N,M)  ,  then  det(^5^5j^  =  0 ,  and  a  pseudo  inverse  does  not  exist.  If  B^is 

full  rank  and  has  a  small  singular  value,  the  inverse  may  be  poorly  eonditioned.  The 
pseudo  inverse  solution  also  known  as  the  Moore-Penrose  inverse  solution  is  used  in 

Equation  (5.6),  where  B^  =B^(^BB^^  ,  whieh  minimizes  the  2-norm  solution  of  the 

following  eonstrained  minimization  problem,  shown  in  Equation  (5.12). 


min ||c||  subjeet  to  B^c  =  5 

c 

II  l|2  T 

where,  c  =  c  c 


(5.12) 


The  pseudo-inverse  is  a  speeial  ease  of  the  weighted  minimum  2-norm  solution 
where  a  weighting  matrix  is  ineluded  in  the  minimization  problem,  resulting  in  the 
following  problem  formulation. 

II  l|2 

min  c  subjeet  to  =  5 

c 

where,  ||c||g  =  Qc_  and  Q  =  Q^  >0 
(5.13) 

B]  =Q-^B^  [bQ-'B^Y  (5.14) 

If  a  matrix  is  not  full  rank  or  if  the  matrix  is  poorly  eonditioned,  an  inverse 
solution  ean  found  using  a  2-norm  and  least  squares  minimization  problem,  where  P  and 
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Q  are  positive  definite  weighting  matriees  (Wie,  2001).  A  singularity  robust  inverse  is 
obtained  using  Equation  (5.16).  For  the  experimental  setup,  the  dimensions  of  P  will  be 
2N  X  2N  and  the  dimensions  of  Q  will  be  M  x  M. 

min|(5^c -5)^  P{B^c_-  s)  +  Qc^  (5.15) 

c  =  B*s_ 

r  n-l  r  n  (5-16) 

where,  B^  =  ^B^PB^  +  Qj  B^P  and  ^B^PB^  +  Qj  is  positive  definite 

Therefore  the  resulting  eontrol  law  ean  be  written  as  (5.17).  In  addition  to 
adjusting  the  gain,  g,  the  weighting  matriees  ean  also  be  adjusted.  The  P  and  Q  matriees 
are  assumed  to  be  diagonal  matriees.  However,  the  diagonal  values  ean  be  adjusted  to 
weight  the  error  and  eontrol.  The  error  is  weighted  using  the  P  matrix  while  the  eontrols 
are  weighted  by  the  Q  matrix.  By  adjusting  the  P  and  Q  matriees  the  stability  of  the 
system  response  ean  be  tuned  experimentally. 

fi,«  =&-s{[b:pB,  +qX  (5.17) 

3.  Direct  Iterative  Zonal  Feedback  Control 

The  iterative  feedbaek  eontrol  law  ean  be  implemented  using  a  direet  zonal 
representation  of  the  wavefront  in  the  influenee  matrix  and  in  the  feedbaek  resulting  in 
Equation  (5.18).  This  method  reduees  the  size  of  the  influenee  matrix  to  a  N  x  P  matrix 
and  reduees  the  eolumn  veetor,  ,  to  N  x  1 .  One  advantage  of  this  method  over  the 

indireet  method  is  that  the  sensor  data  has  a  physieal  meaning  at  eaeh  measured  phase 
point.  This  method  of  wavefront  estimation  is  analogous  to  a  finite  element  approaeh, 
where  if  the  number  of  phase  points  approaehed  infinity  the  wavefront  would  be 
represented  exaetly.  The  ehallenge  with  this  eontrol  algorithm  is  that  the  eontroller  is 
only  as  good  as  the  zonal  phase  estimation.  Typieally  the  A  matrix,  whieh  relates 
measured  slopes  to  phase,  is  very  poorly  eonditioned  resulting  in  a  poor  wavefront 
estimation. 

£,«  =  £,  - (^. )  =  a  - S  (^. )  (5.18) 
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4. 


Direct  Iterative  Modal  Feedback  Control 


The  iterative  modal  feedbaek  eontrol  law  represents  the  influenee  matrix  and 
sensor  data  as  the  modal  eoeffieients  of  the  Zerntke  polynomial  given  as  the  veetor  a^. . 

The  modal  eoeffieients  are  ealeulated  using  the  zonal  phase  representation.  This  method 
eombines  the  indireet  iterative  feedbaek  eontrol  algorithm  with  both  zonal  and  modal 
phase  estimation.  The  addition  of  the  zonal  to  modal  eonversion  adds  additional 
eomputations  to  the  algorithm,  but  also  improves  wavefront  estimation  by  optimally 
fitting  the  zonal  phase  estimation  to  a  polynomial.  One  benefit  to  this  approaeh  is  the 
ability  to  interpret  the  wavefront  easily  as  Zernike  terms.  Additionally  the  influenee 
matrix,  ,  is  often  well  eonditioned  improving  numerieal  stability  when  inverting  the 

influenee  matrix. 

(5.19) 

The  modal  approaeh  ean  be  used  with  a  redueed  number  of  numerieal  ealeulations 
by  applying  the  slope  model  deseribed  in  Equations  (4.17)  to  (4.20).  This  allows  the 
modal  eoeffieients  to  be  eomputed  without  zonal  phase  estimation.  The  slopes  are  related 
to  the  gradient  of  the  Zernike  polynomials  and  the  eoeffieients  are  solved  using  a  least 
squares  approaeh.  This  method  is  more  numerieally  stable  as  there  are  fewer  matrix 
inversions  required  in  the  wavefront  estimation  portion  of  the  eontrol  law. 

Oj^  (5.20) 

C.  ITERATIVE  GRADIENT  FEEDBACK  CONTROL 

The  seeond  eontroller  iteratively  adjusts  the  eontrol  signals  to  reduee  the  varianee 
of  the  wavefront  or  the  varianee  of  the  slope  measurements.  This  is  done  by  taking  the 
derivative  of  the  varianee  with  respeet  to  the  eontrol  signal  to  eompute  an  updated  eontrol 
signal.  This  eontroller  is  similar  to  the  previous  eontroller  exeept  the  gain  is  eomputed  by 
ealeulating  the  derivative  of  the  varianee. 
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Figure  18  Iterative  Gradient  Feedback  Control 

1.  Direct  Gradient  Approach 

A  direct  control  method  developed  by  Zhu,  Sun,  Bartsch,  and  Freeman  begins  by 
representing  the  surface,  ^  of  the  initial  mirror  shape  using  Zernike  polynomials, 

shown  in  Equation  (5.21).  The  initial  surface  is  computed  using  the  Shack-Hartmann 
wavefront  sensor  and  determining  the  Zemtke  coefficients  using  the  wavefront  estimation 
approaches  discussed  in  Chapter  IV.  The  resulting  wavefront  is  described  by  a  vector  of 
Zernike  coefficients.  The  coefficient  vector  a  is  made  up  of  three  other  coefficient 

vectors  as  shown  in  Equation  (5.22).  The  vector  is  the  initial  wavefront  coefficients, 
is  the  input  disturbance  coefficients,  and  is  the  control  coefficients  computed  using 
Equation  (4.16)  or  (4.20),  resulting  in  the  coefficient  vector  in  Equation  (5.24). 

M 

SQ{x,y)  =  Yj(^0k^kix^y)  (5-21) 

k=\ 

a-a^  (5.22) 

^ter=^+^n  (5-23) 

a  =  +  (5.24) 

The  cost  function  of  the  control  algorithm  is  defined  as  the  wavefront  variance 
over  the  entire  aperture.  The  variance  of  the  measured  wavefront  is  related  to  the  control 
signal  so  that  the  variance  can  be  reduced  iteratively,  by  adjusting  the  control  signal.  The 
wavefront  variance  is  described  over  the  unit  circle  for  Cartesian  coordinates  in  Equation 
(5.25).  Assuming  a  planar  reference  wave,  (/)Q{x,y)  ,  the  Zernike  polynomial  expansion 
describes  the  wavefront  variance  from  the  planar  wave,  as  shown  in  Equation  (5.26). 
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Therefore  due  to  the  orthogonality  of  Zemike  polynomials,  when  the  wavefront  varianee 
is  deseribed  by  Equation  (5.26),  the  wavefront  varianee,  cr^,  is  the  sum  of  individual 

varianees  of  the  k*  polynomial  term,  .  The  error  ean  then  be  defined  as  the  varianee, 

as  shown  in  Equation  (5.27). 

I — r  2 

i  (5.25) 


(5.26) 

(5.27) 


The  varianee  equation  ean  be  rewritten  in  terms  of  the  Zemike  polynomial. 
Equation  (5.28).  The  Zemike  eoeffieients  are  eonstants  and  ean  be  moved  in  front  of  the 
integrals.  Sinee  the  Zemike  polynomial  is  valid  over  the  unit  eirele,  the  Zemike  terms 
ean  be  integrated  numerieally  over  the  unit  eirele  resulting  in  a  veetor  of  eoeffieients,  . 
The  error  equation  ean  then  be  rewritten  as  Equation  (5.30),  where  veetors  a  and  w  are 


the  Zemike  eoeffieients  of  the  Zemike  terms  evaluated  over  the  unit  eirele  respeetively. 
Equation  (5.24)  ean  then  be  substituted  into  Equation  (5.30)  resulting  in  the  error 
Equation  (5.31). 


1  Mi—x  1  yji—x 

I  =  j  (^’ t)]  dxdy  =  , 


(5.28) 


where. 

Y  1  Vk?  ^ 

wl=-\  j  [z^{x,y)]dxdy 

(5.29) 

E  =  ^alwl=\a^w\ 

^=1 

(5.30) 

E  =  \{Bc+q^,^^Ywf 

(5.31) 
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Equation  (5.31)  computes  the  error  in  terms  of  a  eontrol  signal,  c,  and  the 
measured  Zemike  eoeffieients.  The  veetor  of  Zernike  terms  is  a  eonstant  that  weights  the 
Zernike  eoeffieients.  The  error  equation  is  a  eost  funetion,  and  the  error  ean  be 
minimized  by  taking  the  derivative  of  the  eost  funetion  with  respeet  to  the  eontrol 
variable,  as  shown  in  Equation  (5.32),  and  the  direetion  of  inereasing  wavefront  error 
with  respeet  to  the  eontrol  signal  ean  be  determined.  Adjusting  the  eontrol  signal  veetor 
in  the  opposite  direetion  of  the  inereasing  wavefront  error  ean  reduee  the  wavefront  error, 
driving  the  gradient  toward  zero.  An  iterative  feedbaek  eontrol  law  ean  be  applied  using 
this  approaeh  by  updating  the  eontrol  voltage  in  Equation  (5.33),  where  //  is  a  sealar 
gain  (Zhu,  Sun,  Bartseh,  Ereeman  &  Eainman,  1999).  The  gain,  // ,  is  determined 
empirieally  and  as  the  eontrol  signals  are  updated  the  error  will  approaeh  a  minimum. 

=  2B^  («  where  =  T wf ,  (5.32) 

dc  ^  '  L  j 

£,+1  (5.33) 

The  above  algorithm  ean  be  optimized  and  tuned  by  adjusting  the  gain,  ju ,  and 

the  veetor  .  Easter  eonvergenee  toward  a  minimum  wavefront  error  ean  be  aehieved 
by  inereasing  the  gain,  ju  ,  however  as  the  algorithm  iterates  a  large  ju  value  may  lead  to 
instabilities.  Therefore  the  gain,  ju ,  ean  be  adaptively  deereased  as  the  error  eonverges 
toward  zero,  improving  stability  by  deereasing  the  step  size  and  preventing  large 
oseillations  around  the  minimum  point.  The  Zernike  term  veetor,  ,  ean  be  adjusted  to 
ehange  the  eonvergenee  speed  of  eaeh  Zernike  term  to  improve  system  performanee. 

One  advantage  to  this  algorithm  is  that  a  modal  disturbanee  ean  be  input  into  the 
system  by  adding  the  term.  This  allows  individual  Zemike  modes  to  be  exeited.  In 

experimental  work,  a  disturbanee  ean  be  easily  added  to  the  system  and  eharaeterized  in 
terms  of  the  Zernike  mode. 
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2. 


Indirect  Gradient  Approach 


An  attractive  control  method  is  to  combine  the  gradient  approach  in  the  previous 
section  with  an  indirect  approach,  eliminating  the  need  to  compute  a  matrix  inverse.  This 
is  accomplished  by  redefining  the  cost  function  used  by  Zhu,  Sun,  Bartsch,  and  Freeman 
(Zhu,  Sun,  Bartsch,  Freeman  &  Fainman,  1999).  Instead  of  defining  error  as  the  variance 
of  the  wavefront  or  the  aperture,  error  is  defined  as  the  sum  of  the  squares  of  individual 
sensor  slopes,  shown  in  Equation  (5.34).  This  becomes  the  cost  function  to  minimize. 

The  slopes  consists  of  both  slopes  in  the  x  and  y  direction,  resulting  in  a  column 
vector  with  2N  rows.  The  slope  vector,  5,  is  defined  by  Equation  (5.35),  where  ^ 

represents  the  initial  slopes  and  ^  the  controlled  slopes.  The  controlled  slopes  can  be 
computed  using  the  influence  matrix  and  control  voltages  as  shown  in  Equation  (5.5). 
The  error  can  be  expressed  as  a  summation  of  the  slopes  in  Equation  (5.36),  and  can  be 
expressed  using  matrices  in  Equation  (5.37),  where  Equation  (5.5)  has  been  substituted 
for  ^  .  By  taking  the  derivative  of  the  error  with  respect  to  the  control  voltage  results  in 

Equation  (5.38),  the  gradient  of  the  slopes  is  determined.  The  derivative  of  the  error  is 
used  in  an  iterative  control  loop  to  determine  the  direction  of  increasing  slope  with 
respect  to  the  control  signal,  as  shown  in  Equation  (5.39),  where  /yis  a  gain  and  ^  is  a 
column  vector  of  zeros  for  a  planar  wave. 


II 

(5.34) 

1  =  ^+^ 

(5.35) 

2N  2 

k=\ 

(5.36) 

(5.37) 

|^  =  2S;(a+S,£)  =  2Sj(a+i) 

(5.38) 
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The  resulting  control  law,  shown  in  Equation  (5.39),  does  not  require  any  inverse 
matrix  operations  or  wavefront  estimation.  However,  a  reference  slope  vector  is  required 
when  computing  the  slopes.  In  this  case  the  reference  slope  vector  is  a  Shack-Hartmann 
sensor  reading  when  the  wavefront  is  planar,  and  therefore  is  a  vector  of  zeros.  The 

algorithm  drives  the  slopes  to  zero  by  changing  the  control  signal,  and  the  resulting 
wavefront  should  match  the  reference  wavefront. 

D,  COMBINED  ITERATIVE  AND  GRADIENT  FEEDBACK 

Another  control  law  is  developed  by  combining  the  direct  iterative  feedback 
control  with  the  gradient  feedback  control,  as  shown  in  Equation  (5.40).  This  control  law 
subtracts  the  wavefront  difference  multiplied  by  a  gain  and  corrects  the  wavefront  by 
compensating  for  the  change  in  the  wavefront  error  by  taking  the  derivative.  This 
algorithm  should  allow  wavefront  error  to  converge  faster,  and  also,  overcome  the 
inherent  local  minimum  problem  with  the  gradient  approach. 

[a^w)  (5.40) 

E.  FILTERING 

Filtering  techniques  can  be  applied  to  control  resonant  peaks  and  known 
frequencies  that  excite  a  structure.  This  is  a  potential  problem  for  a  large  aperture  space 
telescope  as  the  control  system  may  excite  the  structure.  Additionally  a  structure  with 
low  damping  may  lead  to  increased  resonant  peaks  in  the  frequency  response  of  the 
structure  causing  potential  wavefront  error  as  the  aperture  surface  dynamically  changes. 
External  narrowband  disturbances  are  also  of  concern  as  they  can  impart  a  disturbance  on 
a  structure.  Notch  filters  can  cancel  the  narrowband  disturbances.  Notch  filters  are 
characterized  by  a  unit  gain  at  all  frequencies  except  at  the  sinusoidal  frequencies  in 
which  their  gain  is  zero  (Pei  S.  C.  &  Tseng  C.  C.,  1997).’ 

Second  order  structural  filters  can  be  applied  to  compensate  for  known  resonant 

frequencies.  The  second  order  structural  filter  takes  the  form  of  Equation  (5.41),  where 
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the  (o  and  terms  ean  be  manipulated  to  aehieve  phase  lead,  phase  lag,  noteh, 
minimum  bandpass,  non-minimum  phase  all  pass,  and  non-minimum  high  and  low  pass 
filters.  A  seeond  order  noteh  filter  has  the  property  that  all  the  frequeneies  are 


equal,  co^=o)^=co^,  and  the  max  gain  is  determined  by  =  20 log 
2^5 


10 


c 


+  - 


+  1 


G{s)  = 


CO, 


CO, 


r  2^5 


-  + 


+  1 


oo„ 


co„ 


(5.41) 


In  order  to  implement  a  noteh  filter  on  the  experimental  testbed  a  diserete  time 
noteh  filter  is  required.  A  seeond  order  noteh  filter  in  the  Z-domain  is  presented  in 
Equation  (5.42)  where  co^  and  BW  are  the  normalized  eentral  angular  frequeney  and 

bandwidth  of  the  noteh  filter  (Hsu,  Chen  &  Tsai,  2007).  This  noteh  filter  can  be  applied 
to  an  individual  actuator  control  signal,  requiring  a  notch  filter  on  each  control  channel  of 
the  deformable  mirror  to  remove  known  narrowband  disturbances.  The  notch  filter  is 
implemented  on  the  testbed  using  Equation  (5.45)  where  y{n)  represents  the  current 

output  and  u{n)  is  the  current  control  input. 


k2^  Z  +  ^2^ 


(5.42) 


where, 

kj  =-cos(oJ 


1-tan 


K=- 


BW 


1-1- tan 


BW 


(5.43) 

(5.44) 


\  +  k. 


y(n)  =  -k^  (l-i-kj)  j(«-l)-k2j(«-2)-i - -u(^n)  +  k^  (\  +  k2)u{^n-\^ 


(^2+1) 


u{n-2)  (5.45) 
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The  Bode  plot  of  the  noteh  filter  is  shown  in  Figure  19  for  a  frequeney  of  5  Hz 
and  a  bandwidth  of  O.l;^  and  O.S;^.  The  performance  of  the  notch  filter 

on  the  testbed  is  dependent  on  the  a  priori  knowledge  of  the  disturbance  frequency.  The 
disturbance  frequency  input  into  the  system  is  dependent  on  the  sampling  time  of  the 
computer  controller  which  can  vary  in  the  experimental  setup.  Therefore  using  a  narrow 
bandwidth  may  miss  the  uncertain  disturbance  frequency.  Increasing  the  bandwidth  of 
the  notch  filter  will  help  remove  the  disturbance  if  the  frequency  is  not  exactly  known, 
but  will  also  decrease  the  gain  of  the  desired  control  signal. 

The  notch  filter  may  also  be  implemented  in  a  cascade  fashion  to  remove  multiple 
frequencies.  The  Bode  plot  in  Figure  20  shows  a  cascaded  notch  filter  for  0.8  Hz  and  5 
Hz  at  a  bandwidth  of  0.1;^  and  0.5;r.  At  a  higher  bandwidth  of  0.5;r,  the  two  notch 
filters  begin  to  overlap  causing  a  gain  decrease  between  the  two  frequencies.  This  gain 
drop  will  translate  to  a  smaller  control  signal  gain,  reducing  the  overall  control  influence 
of  the  mirror  face-sheet. 


Bode  Diagram 
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Figure  19  Bode  Plot  of  Discrete  Second  Order  Notch  Filter  for  Different  Bandwidths 
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Figure  20  Bode  Plot  of  Cascaded  Discrete  Second  Order  Notch  Filter  for  Different 
Bandwidths 
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VI.  ANALYSIS 


A,  WAVEFRONT  CORRECTION 

The  adaptive  opties  testbed  is  used  to  verily  wavefront  eorreetion  algorithms.  In 
order  to  analyze  the  wavefront  two  figures  of  merit  are  used;  peak  to  valley  wavefront 
aberrations  and  root  mean  square  (RMS)  wavefront  error.  The  RMS  wavefront  error  and 
peak  to  valley  wavefront  aberration  is  easily  eomputed  from  a  modal  phase 
representation.  Therefore  to  eompare  results  of  different  eontrol  algorithms  a  standard 
wavefront  estimation  method  is  used  when  eomputing  error  and  reeonstrueting  the 
wavefront. 

The  standard  wavefront  estimation  algorithm  used  is  the  modal  phase  estimation 
method  deseribed  in  Equation  (4.20)  using  the  derivatives  of  the  Zernike  polynomial 
terms.  This  method  was  ehosen  beeause  the  wavefront  ean  be  direetly  estimated  from  the 
Shaek-Hartmann  slope  data.  The  wavefront  estimation  method  only  requires  one  matrix 
inverse  of  a  well  eonditioned  matrix  of  orthogonal  Zernike  terms.  A  total  of  21  Zernike 
terms  are  used  in  the  wavefront  estimation  to  eompute  the  RMS  wavefront  error  and  peak 
to  valley  wavefront  aberration. 

1,  Peak  to  Valley  Wavefront  Aberation 

When  eomplieated  wavefront  aberrations  are  deseribed  in  terms  of  third  order 
aberrations,  the  piston,  tip,  and  tilt  aberration  is  ignored.  The  easiest  way  to  deseribe  the 
wavefront  error  is  in  terms  of  peak  to  valley  wavefront  aberrations,  measured  in  waves 
(Wyant,  2003).  This  figure  of  merit  deseribes  how  far  the  measured  wavefront  is  from 
the  desired  wavefront.  This  figure  of  merit  is  an  absolute  value  of  the  measured 
differenee  from  the  desired  wavefront,  similar  to  an  infinity  norm.  The  peak  to  valley 
wavefront  aberration  is  ealeulated  using  a  wavefront  estimation  teehnique  and  summing 
the  absolute  value  of  the  maximum  and  minimum  deviation  from  a  planar  wavefront. 
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2. 


RMS  Wavefront  Error 


The  use  of  the  peak  to  valley  figure  of  merit  ean  be  misleading  as  it  does  not 
deseribe  the  area  over  whieh  the  error  oeeurs  (Wyant,  2003).  The  peak  to  valley 
measurement  only  gives  the  maximum  wavefront  error.  Another  figure  of  merit  that 
provides  a  more  meaningful  measure  of  wavefront  quality  is  the  root  mean  square  (RMS) 
wavefront  error.  The  RMS  wavefront  error  is  the  square  root  of  the  varianee  of  the 
wavefront.  The  varianee  was  deseribed  in  ehapter  V  and  used  as  a  eost  funetion. 
Therefore  by  taking  the  square  root  of  the  double  integral  of  the  Zernike  polynomial  over 
the  unit  eirele  the  RMS  error  is  eomputed,  using  Equation  (6.1). 


cr^  = 


1  vi- 


n 


J  J  [a,,z^[x,y)\dxdy  = 


j  1  WA  _ _ 

I  " V 


(6.1) 


B,  SURFACE  CORRECTION  OF  THE  MMDM 

The  eontrol  algorithms  developed  previously  were  applied  to  the  MMDM  to 
eorreet  the  mirror  surfaee  from  a  biased  position.  This  allows  the  eontrol  algorithms  to 
be  validated  and  eompared  when  not  subjeet  to  a  disturbanee.  A  bias  eontrol  signal  of 
240  was  applied  to  all  the  aetuators  on  the  MMDM  resulting  in  a  peak  to  valley 
wavefront  aberration  of  17.7  waves  and  RMS  wavefront  error  of  10.96  as  shown  in 
Figure  21.  In  all  experimental  tests,  the  eontrol  signals  were  limited  between  0  and  240 
based  on  the  results  from  Figure  15.  The  error  of  the  wavefront  is  only  oaleulated  over 
85%  of  the  Hartmann  mask  aperture  using  a  diameter  of  3  mm.  This  eliminates  the  outer 
fringes  of  the  wavefront  as  the  MMDM  deformation  range  is  limited  at  the  outer  edges 
sinee  the  membrane  is  fixed  at  the  boundary.  The  eontrol  algorithms  were  eonfigured  to 
drive  the  mirror  from  the  biased  position  to  a  position  resulting  in  a  planar  wavefront. 
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Figure  21  Biased  Wavefront  of  MMDM,  cr  =  10.96  ,  Peak  to  Valley  =  17.79  waves, 
Control  Signal  Equals  240 

1,  Indirect  Iterative  Control 

The  indirect  iterative  control  methods  include  the  control  algorithms  presented  in 
Equations  (5.11)  and  (5.17).  The  gain  is  set  to  0.9  in  both  experiments  and  the  loop  is  run 
for  150  iterations.  Eigure  22  and  Eigure  23  show  the  wavefront  representation  of  the 
Ecorrected  wavefront  after  150  iterations.  Both  methods  improve  the  biased  wavefront, 
however  the  singularity  robust  method  outperforms  the  SVD  method  by  a  factor  of  10. 
Both  control  algorithms  are  the  same  except  for  the  matrix  inversion  algorithm.  Eigure 
24  shows  the  RMS  wavefront  error  history  for  the  150  iterations  for  both  algorithms. 

The  numerical  instabilities  of  inverting  the  254  x  37  influence  matrix  using  the 
SVD  inverse  are  apparent  in  the  error  history  plot.  The  jagged  peaks  seen  in  the  error 
using  the  SVD  method  in  Eigure  24  indicate  an  over  correction  and  instability.  This  can 
be  improved  by  decreasing  the  gain.  The  singularity  robust  method  uses  the  weighting 
matrices  Q  and  P  to  weight  the  influence  of  the  control  and  the  error.  In  this  experiment 
the  Q  matrix  is  a  diagonal  matrix  of  ones  with  the  size  37  x  37,  and  the  P  matrix  is  a 
diagonal  matrix  of  0.01  with  the  size  of  254  x  254.  The  stability  and  performance  of  the 
control  loop  are  improved  by  using  the  singularity  robust  inverse. 
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Figure  22  Corrected  Wavefront  Using  Indirect  Iterative  Feedback  Control  with 
Shack-Hartmann  Slopes  and  SVD  Inverse,  Peak  to  Valley  =  5.42  Waves 


0.2 


Figure  23  Corrected  Wavefront  Using  Indirect  Iterative  Feedback  Control  with 
Shack-Hartmann  Slopes  and  Singularity  Robust  Inverse,  Peak  to  Valley  = 
0.266  Waves 
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Figure  24  Error  Flistory  Using  Indirect  Iterative  Feedback  Control  with  Shack- 
Hartmann  Slopes  with  Singularity  Robust  Inverse  and  SVD  Inverse 

2,  Direct  Iterative  Control 

Three  direct  iterative  control  methods  are  employed  with  the  difference  being  the 
phase  representation.  The  control  algorithm  in  (5.18)  uses  a  zonal  phase  representation 
where  the  wavefront  is  computed  from  the  Shack-Hartmann  slope  measurements  using  a 
bi-quadratic  spline  producing  the  relationship  in  Equations  (4.10)  and  (4.11)  which  are 
modified  from  the  Southwell  geometry  in  Equations  (4.8)  and  (4.9)  (Southwell,  1980). 
The  resulting  wavefront  and  error  history  are  shown  in  Figure  25  and  Figure  26.  The 
zonal  iterative  feedback  control  algorithm  does  show  improvement  from  the  biased 
wavefront,  but  still  does  a  poor  job.  This  is  due  to  the  poor  condition  number, 
K  w  1x1 0'^ ,  of  the  A  matrix  that  computes  the  phase  difference  between  two  grid  points  in 
Equation  (4.13).  Therefore,  the  feedback  used  in  the  control  algorithm,  (j) ,  is  not 

numerically  stable  and  accumulates  error.  The  gain  for  the  control  loop  is  set  to  0.2 
empirically  to  provide  a  stable  and  fast  response. 
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Wavefront  Surface,  Peak  to  Valley  =  6.0926 
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Figure  25  Corrected  Wavefront  Using  Zonal  Iterative  Feedback  Control,  Peak  to 
Valley  =  6.0926  Waves 


Figure  26  Error  Flistory  Using  Zonal  Iterative  Feedback  Control 

The  control  algorithm  in  Equation  (5.19)  uses  modal  wavefront  estimation 
derived  from  zonal  wavefront  estimation.  The  Zemike  coefficients  are  computed  using 
(4.16)  and  are  used  as  the  feedback  to  the  control  algorithm.  This  method  performs 
significantly  better  than  the  previous  zonal  method.  The  error  present  in  the  zonal 
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wavefront  estimation  is  redueed  as  the  zonal  wavefront  is  fitted  by  a  least  squares 
teehnique  to  21  Zemike  eoeffieients.  The  gain  for  this  method  was  set  to  0.4  empirically. 
This  gain  produced  good  results  with  a  stable  response  and  fast  convergence. 


0.2 


y.  mm 


Figure  27  Corrected  Wavefront  Using  Modal  (derived  from  Zonal)  Iterative 
Feedback  Control,  Peak  to  Valley  =  0.268  Waves. 


Figure  28  Error  History  Using  Modal  (derived  from  Zonal)  Iterative  Feedback 
Control 
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The  third  direct  iterative  control  method  estimates  the  wavefront  directly  from  the 
Shack-Hartmann  wavefront  sensor  slope  measurements  using  the  derivatives  of  the 
Zernike  terms.  The  control  algorithm  in  Equation  (5.20)  outperforms  the  two  previous 
control  methods.  The  gain  was  set  to  0.11  empirically.  This  provided  a  stable  response 
that  converged  relatively  fast.  The  improved  performance  comes  from  the  simplified 
wavefront  estimation  technique.  The  least  squares  wavefront  estimation  requires  the 
matrix  inversion  of  the  matrix  dZ ,  which  is  well  conditioned  since  it  contains  the 
orthogonal  Zemike  terms,  reducing  numerical  error.  The  resulting  influence  matrix  is 
also  well  conditioned  making  the  control  loop  numerically  stable. 


Figure  29  Corrected  Wavefront  Using  Modal  (derived  from  Zemike  Derivatives) 
Iterative  Feedback  Control,  Peak  to  Valley  =  0.0826 
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Figure  30  Error  History  Using  Modal  (derived  Zemike  Derivatives)  Iterative 
Feedback  Control 

3,  Iterative  Gradient  Feedback  Control 

Two  iterative  gradient  feedback  control  algorithms  were  implemented  using 
Equations  (5.33)  and  (5.39).  The  first  control  algorithm  iteratively  adjusts  the  control 
voltage  by  minimizing  the  variance  of  the  wavefront,  estimated  using  the  modal 
estimation  method  described  in  (4.20)  (Zhu,  Sun,  Bartsch,  Freeman  &  Fainman,  1999). 
The  gain  is  empirically  set  to  0.039.  The  results  in  Figure  31  and  Figure  32  show  a 
wavefront  improvement  over  the  150  iterations. 

The  other  iterative  gradient  feedback  control  algorithm  was  developed  by 
minimizing  the  measured  slope  error.  This  algorithm  does  not  require  wavefront 
estimation  and  can  be  implemented  directly  from  the  Shack-Hartmann  wavefront  sensor 
feedback.  The  gain  is  empirically  set  to  0.00001.  The  results  in  Figure  33  and  Figure  34 
show  wavefront  improvement  over  the  150  iterations  outperforming  the  variance 
minimization  method. 
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Figure  31  Corrected  Wavefront  Using  Iterative  Gradient  Feedback  Minimizing 
Wavefront  Variance,  Peak  to  Valley  =  0.735  Waves 


Figure  32  Error  History  Using  Iterative  Gradient  Feedback  Minimizing  Wavefront 
Variance 
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Figure  33  Corrected  Wavefront  Using  Iterative  Gradient  Feedback  Minimizing 
Slope  Error,  Peak  to  Valley  =  0.33714 


Figure  34  Error  History  Using  Iterative  Gradient  Eeedback  Minimizing  Slope  Error, 
150  Iterations 

The  iterative  gradient  feedback  control  algorithms  performance  can  be  improved 
by  adjusting  the  gain  and  adjusting  the  Zemtke  coefficient  weighting  vector,  ,  for 
Equation  (5.33).  This  can  improve  the  convergence  speed  as  well  as  the  wavefront  error 


55 


in  both  algorithms.  One  inherent  problem  with  gradient  methods  is  that  a  global 
minimum  is  not  guaranteed.  Therefore,  the  algorithm  may  drive  the  wavefront  to  a  loeal 
minimum,  limiting  the  performanee  of  the  algorithm. 

In  the  two  proposed  gradient  methods  the  gain  must  be  determined  empirieally  to 
eontrol  the  step  size  of  eontrol  voltage  ehange  for  eaeh  iteration.  If  the  step  size  is  too 
large  an  oseillatory  response  may  develop  as  the  mirror  oseillates  around  a  loeal 
minimum.  Using  the  previous  iterative  feedbaek  methods,  the  gain  is  between  0  and  1 
beeause  the  eontrol  law  is  based  on  an  influenee  matrix  that  is  built  using  a  eontrol  signal 
between  0  and  1,  making  the  eontroller  easier  to  tune. 

4,  Combined  Iterative  and  Gradient  Feedback  Control 

The  direet  iterative  feedbaek  eontrol  law  using  Zernike  modes  estimated  from 
Zernike  derivatives  is  eombined  with  the  direet  gradient  feedbaek  minimizing  the 
varianee  eontrol  law.  The  results  were  eompared  against  the  direet  iterative  feedbaek 
eontrol  law  using  Zernike  eoeffieients  estimated  from  Zernike  derivatives.  The  gain  for 
the  direet  iterative  feedbaek  was  set  to  0.4  and  the  gain  for  the  gradient  feedbaek  was  set 
to  0.01.  The  eombined  eontrol  law  outperformed  the  standard  iterative  modal  feedbaek 
eontrol  law  as  seen  in  Figure  35,  Figure  36,  and  Figure  37.  The  RMS  wavefront  error 
and  peak  to  valley  wavefront  error  is  over  three  times  smaller.  The  wavefront  error  also 
eonverges  faster. 
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Figure  35  Corrected  Wavefront  Using  Modal  (derived  from  Zemike  Derivatives) 
Iterative  Feedback  Control,  Peak  to  Valley  =  0.0788 
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Figure  36  Corrected  Wavefront  Using  Modal  (derived  from  Zernike  Derivatives) 
Iterative  Feedback  with  Gradient  Feedback  Control,  Multiple  Gain,  Peak 
to  Valley  =  0.0219 
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Figure  37  Error  History  Comparing  Iterative  Modal  Feedback  and  Iterative  Modal 
Feedback  Combined  with  Gradient  Feedback 

5,  Control  Method  Comparison 

The  wavefront  estimation  process  of  wavefront  control  makes  up  a  large  portion 
of  the  control  loop.  Although  modal  phase  estimation  adds  an  additional  process  to  the 
wavefront  control  algorithm  over  an  indirect  method,  it  can  also  decrease  the  overall 
number  of  computations  required  per  iteration.  The  number  of  numerical  steps  using 
modal  phase  estimation  is  dependent  on  the  number  of  Zernike  terms  used  to  describe  the 
wavefront.  Therefore,  it  is  possible  to  reduce  the  number  of  numerical  computations  in 
the  control  algorithm  by  using  a  modal  phase  estimation  technique  rather  than  using  an 
indirect  wavefront  representation.  Figure  38  shows  the  number  of  multiplication 
operations  required  for  one  iteration  vs.  the  number  of  Zernike  terms.  The  indirect 
method  remains  constant  as  it  is  dependent  on  the  number  of  lenslets  in  the  Hartmann 
mask,  while  the  direct  modal  method  is  linearly  dependent  on  the  number  of  Zernike 
terms.  In  the  previous  experiments  only  21  Zernike  terms  are  used,  therefore  the  direct 
modal  iterative  control  algorithm  has  fewer  multiplications  than  the  indirect  method 


58 


using  slope  information.  However,  if  the  number  of  Zemike  terms  is  inereased  past  34  it 
is  eomputationally  benefieial  to  use  a  control  algorithm  with  an  indirect  wavefront 
representation. 


Figure  38  Multiplication  Operations  For  Single  Control  Loop  Iteration  vs.  Number 
of  Zemike  Modes  Used,  Using  37  Actuator  Mirror,  and  127  Lenslet 
Hartmann  Mask,  Modal  Phase  Representation  from  Zemike  Derivatives, 
and  Slope  Representation  from  Measured  x  and  y  Slopes 

Table  5  shows  the  error  for  the  previously  discussed  control  techniques. 
Analyzing  the  non-gradient  iterative  feedback  control  algorithms  first,  shows  that  the 
direct  iterative  feedback  control  using  modal  wavefront  estimation  from  the  Zemike 
derivatives  performed  the  wavefront  correction  the  best.  The  indirect  iterative  feedback 
controller  using  a  singularity  robust  inverse  also  performed  remarkably  well.  Both  these 
techniques  minimized  the  steps  required  to  estimate  the  wavefront,  or  neglected  a 
wavefront  estimation  altogether.  Both  techniques  are  singularity  robust  providing 
stability  when  inverting  the  influence  matrix  and  computing  the  updated  control  signal. 
This  accounts  for  their  superior  performance  during  the  experiment. 

Comparing  the  gradient  approach,  the  slope  minimization  performed  better  than 
the  variance  minimization  method.  This  is  due  to  the  variance  minimization  taking  place 
over  the  continuous  aperture  using  wavefront  estimation  while  the  slope  minimization 
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process  only  takes  plaee  at  the  zonal  phase  points  determined  by  the  Hartmann  mask. 
The  slope  minimization  algorithm  is  likened  to  a  finite  element  model,  where  the 
performanee  could  be  potentially  improved  by  adding  additional  elements  or  lenslets  over 
the  aperture. 

The  best  performanee  is  aehieved  by  eombining  the  direet  iterative  feedbaek 
using  modal  phase  estimation  from  Zernike  derivatives  with  the  gradient  feedbaek 
method  minimizing  the  varianee.  This  method  overeomes  the  inherent  loeal  minimum 
problem  with  gradient  minimization  approaehes  and  also  allows  faster  eonvergenee  by 
adaptively  adjusting  the  eontrol  signal  based  on  the  rate  at  whieh  the  error  is  ehanging. 


Control  Algorithm 

Peak  to  Valley 

RMS  Wavefront  Error 

Indirect  Iterative  Feedback  SVD 

5.424 

2.819 

Indirect  Iterative  Feedback  SR 

0.266 

0.143 

Direct  Iterative  Feedback,  Zonal 

6.092 

2.383 

Direct  Iterative  Feedback,  Modal  From  Zonal 

0.268 

0.086 

Direct  Iterative  Feedback,  Modal  from  Zernike  Derivatives 

0.082 

0.028 

Iterative  Gradient  Feedback,  Variance  Minimization 

0.736 

0.199 

Iterative  Gradient  Feedback,  Slope  Minimization 

0.337 

0.089 

Combined  Direct  Iterative  and  Gradient  Feedback 

0.022 

0.008 

Table  5  Error  Comparison  for  Wavefront  Control  Methods 


C.  MMDM  CORRECTION  SUBJECT  TO  A  DISTURBANCE 

The  next  experiment  subjeets  the  MMDM  to  a  low  frequeney  sinusoidal 
disturbanee.  The  disturbanee  is  input  into  the  mirror  by  adjusting  the  aetuator  eontrol 
voltages.  This  experiment  is  eondueted  using  the  MMDM.  A  disturbanee  signal  with  a 
frequency  of  5  Hz  and  amplitude  of  54  was  introdueed  to  the  mirror.  The  mirror  is  biased 
at  240  for  the  start  of  eaeh  experiment  and  the  eontroller  is  operated  for  250  iterations. 
Five  of  the  eontrol  laws  developed  in  the  previous  seetion  were  evaluated  for  their  ability 
to  eorreet  a  dynamie  disturbanee  in  the  mirror. 

The  five  eontrol  laws  evaluated  are  the  Indireet  Iterative  Feedbaek  using  a 

singularity  robust  inverse,  direet  iterative  feedbaek  using  modal  phase  estimation  from 

Zernike  derivatives,  the  iterative  gradient  feedbaek  using  varianee  minimization,  the 

iterative  gradient  feedbaek  using  slope  minimization,  and  the  eombined  iterative  and 
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gradient  feedback.  Figure  40  through  Figure  44  show  the  RMS  wavefront  error  history 
for  the  five  controllers  with  and  without  a  disturbance.  Figure  39  shows  the  error 
imparted  on  a  planar  wave  when  the  mirror  is  subjected  to  the  disturbance  signal  with  no 
control  resulting  in  a  peak  RMS  wavefront  error  of  0.8546.  The  peak  steady  state  RMS 
wavefront  errors  are  used  to  compare  the  algorithms,  Table  6.  The  gradient  methods 
performed  poorly  adding  additional  error  to  the  wavefront.  The  iterative  feedback 
methods  including  the  combined  iterative  feedback  with  gradient  feedback  method  were 
all  comparative  in  performance,  driving  down  the  wavefront  error  from  the  biased 
position  but  leaving  residual  error  due  to  the  disturbance  signal.  All  the  controllers  could 
benefit  from  using  notch  filters  to  remove  the  known  sinusoidal  disturbance. 

The  poor  performance  of  the  gradient  method  is  due  to  the  bandwidth  of  the 
sensors  and  the  control  loop,  approximately  30Hz.  The  control  signal  is  computed  based 
on  the  derivative  of  the  error  with  respect  to  the  control  signal.  This  causes  the  controller 
to  over  correct  adding  additional  error  to  the  wavefront.  Increasing  the  bandwidth  of  the 
sensor  and  the  control  loop  would  improve  the  iterative  gradient  controller  performance, 
by  decreasing  the  iterative  step  size.  The  additional  error  may  also  be  reduced  by 
decreasing  the  gain  at  the  cost  of  increasing  the  settling  time. 


Figure  39  RMS  Wavefront  Error  of  Planar  Surface  Subjected  to  5  Hz  Sinusoidal 
Disturbance  with  Amplitude  of  54  Volts 
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Figure  40  Error  History  Using  Indirect  Iterative  Feedback  with  Singularity  Robust 
Inverse,  with  and  without  5  Hz  Sinusoidal  Disturbance 


Figure  41  Error  History  Using  Direct  Iterative  Feedback  with  Modal  Phase 
Estimation  from  Zemike  Derivatives,  with  and  without  5  Hz  Sinusoidal 
Disturbance 
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Figure  42  Error  History  Using  Iterative  Gradient  Feedbaek  Minimizing  Varianee, 
with  and  without  5  Hz  Disturbance 


Figure  43  Error  History  Using  Iterative  Gradient  Feedback  Minimizing  Slopes,  with 
and  without  5  Hz  Disturbance 
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Figure  44  Error  History  Using  Combined  Direct  Iterative  and  Gradient  Feedback, 
with  and  without  5  Hz  Sinusoidal  Disturbance 


No  Disturbance 

Disturbance 

Control  Method 

a 

Umax  at  steady  state 

Indirect  Iterative  Feedback  SR 

0.0176 

0.91261 

Direct  Iterative  Feedback,  Modal  from  Zemike  Derivatives 

0.2317 

0.8897 

Iterative  Gradient  Feedback,  Variance  Minimization 

0.10174 

1.95 

Iterative  Gradient  Feedback,  Slope  Minimization 

0.0288 

1.6483 

Combined  Direct  Iterative  and  Gradient  Feedback 

0.00619 

0.8547 

Table  6  Comparison  of  Control  Methods  Subject  to  a  5  Hz  Sinusoidal  Disturbance 


D.  NOTCH  FILTERING 

1,  Filtering  Noisy  Actuators 

The  discrete  time  notch  filter  is  tested  by  introducing  a  sinusoidal  disturbance  to 
an  individual  actuator  on  the  MMDM  and  applying  the  notch  filter  to  the  control  law  to 
remove  the  disturbance.  A  single  actuator  is  tested  first  to  minimize  coupling  between 
notch  filters  on  multiple  actuators.  The  notch  filter  was  tested  using  the  indirect  iterative 
feedback  control  algorithm  using  a  singularity  robust  inverse.  This  control  algorithm  was 
chosen  to  test  the  notch  filter  because  the  algorithm  runs  consistently  faster  than  the 
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algorithms  with  modal  phase  estimation  and  performs  eomparably  well.  The  control 
algorithm  operates  between  28  and  30  Hz.  With  this  sampling  rate  the  notch  filter  should 
be  able  to  handle  a  low  frequency  sinusoidal  disturbance  less  than  5  Hz. 

a.  5  Hz  Sinusoidal  Disturbance  on  Single  Actuator 

The  notch  filter  was  tested  against  a  5  Hz  sinusoidal  disturbance  with 
amplitude  of  54.  The  notch  filter  was  tuned  to  5  Hz  and  a  bandwidth  of  O.Sti.  Figure  45 
shows  the  RMS  error  history  for  the  surface  of  the  wavefront  with  and  without  the  notch 
filter.  The  notch  filter  reduced  the  sinusoidal  disturbance  and  the  maximum  RMS  error. 
Figure  46  shows  the  control  history  showing  that  the  error  signal  is  still  present  in  the 
control  signal,  but  at  a  significantly  lower  amplitude.  The  remaining  error  is  due  to 
uncertainties  in  the  sampling  rate  changing  the  disturbance  frequency.  Increasing  the 
bandwidth  of  the  notch  filter  will  improve  performance  by  reducing  the  sinusoidal 
disturbance  on  the  actuator,  but  will  also  reduce  the  control  influence  of  the  actuator. 
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Figure  45  Error  History  for  5  Hz  Sinusoidal  Disturbance  with  an  Amplitude  54 
Volts,  on  MMDM  Actuator  10,  with  and  without  Notch  Filter,  Initial  Bias 
of  127 
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Figure  46  Control  History  for  MMDM  Actuator  10  with  and  without  Notch  Filter,  5 
Hz  Sinusoidal  Disturbance,  54  Amplitude 


b.  2  Hz  Sinusodial  Disturbance  on  Single  Actuator 

Using  the  same  actuator  but  applying  a  2  Hz  sinusoidal  disturbance  the 
bandwidth  of  the  notch  filter  is  set  to  O.Iti  and  0.5n.  The  error  history.  Figure  47,  shows 
that  the  RMS  wavefront  error  is  improved  by  setting  the  notch  filter  to  a  higher 
bandwidth.  The  control  history  also  shows  a  reduced  sinusoidal  disturbance  in  the 
control  signal  at  the  higher  bandwidth.  Figure  48.  The  wider  bandwidth  shows 
improvement  because  it  also  reduces  the  gain  of  adjacent  frequencies,  which  helps 
account  for  the  uncertainty  in  the  sampling  time  of  the  computer  control  loop  which  adds 
uncertainty  to  the  disturbance  signal. 
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Figure  47 


Figure  48 
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Error  History  for  2  Hz  Sinusoidal  Disturbanee  with  an  Amplitude  54 
Volts,  on  MMDM  Aetuator  10,  with  and  without  Notch  Filter,  Initial  Bias 
of  127 
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Control  History  for  MMDM  Channel  10  with  and  without  Notch  Filter,  2 
Hz  Sinusoidal  Disturbance,  54  Amplitude 
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c.  Cascaded  Notch  Filters  on  Single  Actuator 

Two  disturbances  at  different  frequeneies  were  introduced  into  actuator  10 
and  a  cascading  notch  filter  was  applied  to  remove  the  disturbanee.  The  easeading  noteh 
filters  were  set  to  0.8  Hz  and  2  Hz  respeetively  with  a  O.Sti  bandwidth  and  suecessfully 
reduced  the  RMS  wavefront  error.  A  residual  disturbance  is  still  evident  in  both  the  error 
and  eontrol  history  seen  in  Figure  49  and  Figure  50.  The  overall  RMS  wavefront  error  is 
greater  than  the  previous  experiments  where  only  one  disturbance  frequeney  was 
introduced.  This  is  partially  due  to  the  overall  gain  drop  at  adjacent  frequeneies  as  seen 
in  the  Bode  plot,  Figure  20. 
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Figure  49  Error  History  for  .08  and  2  Hz  Sinusoidal  Disturbance  with  an  Amplitude 
54  Volts,  on  MMDM  Aetuator  10,  with  and  without  Caseading  Notch 
Filter,  Initial  Bias  of  127 
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Figure  50  Control  Flistory  for  MMDM  Channel  10  with  and  without  Cascading 
Notch  Filter,  0.8  and  2  Hz  Sinusoidal  Disturbance,  54  Amplitude 


d.  2  Hz  and  0.8  Hz  Sinusoidal  Disturbance  on  Non-Adjacent 
Actuators 

A  second  disturbance  was  added  to  the  mirror  at  actuator  16.  Actuators  10 
and  16  are  on  opposite  sides  of  the  face  sheet  to  minimize  coupling.  Both  actuators  were 
subjected  to  a  sinusoidal  disturbance  with  two  different  frequencies,  2  Hz  on  actuator  10 
and  0.8  Hz  on  actuator  16.  The  notch  filters  were  tuned  to  2.0  and  0.8  Hz  respectfully 
and  a  bandwidth  of  0.57i.  The  notch  filters  performed  well  reducing  the  overall  RMS 
wavefront  error,  and  removing  the  sinusoidal  disturbance,  as  shown  in  Figure  51.  The 
amplitude  of  the  disturbance  is  reduced  in  both  control  signals  with  the  higher  frequency 
disturbance  being  reduced  more  than  the  lower  frequency  disturbance,  as  shown  in  Figure 
52. 
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Figure  51  Error  Flistory  for  2.0  FIz  and  0.8  FIz  Sinusoidal  Disturbance  with  an 
Amplitude  54  Volts,  on  MMDM  Actuator  10  and  16,  with  and  without 
Notch  Filter,  Initial  Bias  of  127 


Figure  52  Control  Flistory  for  MMDM  Channel  10  and  16  with  and  without  Notch 
Filter,  2  Hz  Sinusoidal  Disturbance  on  Channel  10  and  0.8  Hz  on  Channel 
16  with  an  Amplitude  54 
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e.  2  Hz  and  0. 8  Hz  Sinusoidal  Disturbance  on  Adjacent  Actuators 

To  explore  the  eoupling  between  aetuator  influenees  on  the  face  sheet  two 
disturbances  were  introduced  to  two  adjacent  actuators,  10  and  11.  The  same  frequency 
disturbance  is  introduced  on  both  actuators  but  with  only  one  notch  filter  on  actuator  10. 
The  overall  RMS  wavefront  error  was  decreased;  however,  a  residual  error  is  still  present 
at  2  Hz  due  to  the  error  caused  by  actuator  1 1,  as  shown  in  Figure  53.  The  control  signal 
history  of  actuator  10  shows  that  the  disturbance  was  successfully  decreased  while  the 
control  history  for  actuator  11  shows  a  decrease  in  the  maximum  control  voltage,  but 
similar  absolute  amplitude  at  steady  state.  The  change  in  control  voltage  demonstrates 
the  coupling  of  the  actuators  due  to  the  continuous  surface  of  the  face  sheet. 
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Figure  53  Error  History  for  2.0  Hz  Sinusoidal  Disturbance  with  an  Amplitude  54 
Volts,  on  MMDM  Actuator  10  and  11,  with  and  without  Notch  Filter, 
Initial  Bias  of  127 
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Figure  54  Control  Flistory  for  MMDM  Channel  10  and  11  with  and  without  Noteh 
Filter,  2  Flz  Sinusoidal  Disturbance  on  Channel  10  and  Channel  1 1  with  an 
Amplitude  54 

A  similar  experiment  to  the  previous  is  conducted  using  two  different 
frequency  disturbances  on  the  two  separate  actuators,  10  and  11,  with  a  notch  filter  only 
on  actuator  10.  The  results  are  similar  to  the  previous  results  with  the  higher  frequency 
disturbance  being  removed  from  actuator  10  and  the  lower  frequency  disturbance  present 
in  both  the  error  and  control  history  plots  due  to  actuator  11,  as  shown  in  Figure  55  and 
Figure  56.  Again,  the  adjacent  notch  filter  impacts  the  control  signal  of  actuator  1 1  when 
implemented  with  the  control  law. 
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Figure  55  Error  History  for  2.0  Hz  and  0.8  Hz  Sinusoidal  Disturbance  with  an 
Amplitude  54  Volts,  on  MMDM  Actuator  10  and  11,  with  and  without 
Notch  Filter  on  Actuator  10,  Initial  Bias  of  127 
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Figure  56  Control  History  for  MMDM  Channel  10  and  11  with  and  without  Notch 
Filter,  2  Hz  and  0.8  Hz  Sinusoidal  Disturbance  on  Channel  10  and 
Channel  1 1  with  an  Amplitude  54 
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Two  notch  filters  were  applied  to  the  adjacent  actuators,  10  and  11.  The 
notch  fdters  were  tuned  to  the  appropriate  frequencies  with  a  bandwidth  of  O.Sti.  The 
results  show  that  the  disturbance  is  notched  out  successfully  reducing  the  RMS  wavefront 
error.  The  error  is  reduced  by  a  factor  of  four  over  using  just  the  iterative  control  law 
alone.  The  improved  RMS  wavefront  error  over  the  previous  experiment  where  a  non- 
adjacent  actuator  was  used  is  due  to  the  localization  of  the  disturbance  on  the  face  sheet. 
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Figure  57  Error  Flistory  for  2.0  Flz  and  0.8  Flz  Sinusoidal  Disturbance  with  an 
Amplitude  54  Volts,  on  MMDM  Actuator  10  and  11,  with  and  without 
Notch  Filter,  Initial  Bias  of  127 
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Figure  58  Control  Flistory  for  MMDM  Channel  10  and  11  with  and  without  Noteh 
Filter,  2  Hz  Sinusoidal  Disturbance  on  Channel  10  and  0.85  Hz  on 
Channel  1 1  with  an  Amplitude  54 


f.  2  Hz  Sinusoidal  Disturbance  on  37  MMDM  Actuators 

In  this  experiment,  a  2  Hz  sinusoidal  disturbance  and  a  discrete  time  notch 
filter  are  applied  to  all  the  actuators.  The  filter  bandwidth  is  set  to  O.Iti.  The  notch  filter 
performed  well  reducing  the  maximum  steady  state  wavefront  disturbance  by  a  factor  of 
four.  The  narrow  bandwidth  of  the  notch  filter  allows  for  a  quick  settling  time  but  allows 
some  of  the  disturbance  to  pass  due  to  the  uncertainty  of  the  disturbance  frequency. 
Increasing  the  bandwidth  will  decrease  the  overall  error  but  will  also  increase  the  settling 
time  by  reducing  the  overall  actuator  gains. 


75 


7 


D 

m 

5 

3 


4  5 

time  (sec) 


Disturbance 
Disturbance  with  NF 


Figure  59  Error  History  for  2.0  Hz  Sinusoidal  Disturbance  with  an  Amplitude  54 
Volts,  on  37  MMDM  Actuators,  with  and  without  Notch  Filter,  Initial  Bias 
of  127 


2,  Filtering  a  Simulated  Vibration 

A  simulated  disturbance  due  to  vibration  is  introduced  into  the  MMDM  using  the 
previously  discussed  modal  iterative  feedback  control  law  using  wavefront  estimation 
from  Zemike  derivatives.  A  disturbance  is  introduced  into  the  mirror  by  introducing  a 
sinusoidal  disturbance  into  the  third  Zemike  term,  which  represents  focus.  This  creates 
an  oscillating  focus  aberration  on  the  surface  of  the  mirror  and  in  the  wavefront.  The 
disturbance  simulates  a  first  mode  vibration  disturbance  across  the  face  sheet  at  a  specific 
frequency.  Unlike  the  previous  notch  filter  experiment  where  each  actuator  had  the  same 
disturbance  signal,  the  actuator  voltages  are  adjusted  to  achieve  the  desired  wavefront 
and  mirror  aberration.  The  second  order  notch  filtering  techniques  previously  developed 
are  applied  to  the  modal  iterative  feedback  control  law  to  reduce  the  disturbance. 

The  resulting  error  history  for  the  disturbance  shows  that  the  notch  filter  was  able 
to  remove  a  portion  of  the  disturbance  successfully  using  a  bandwidth  of  0.27i,  Figure  60. 
However  by  increasing  the  bandwidth  to  O.Sti  the  notch  filter  was  able  to  almost 
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completely  remove  the  simulated  structural  disturbance  across  the  face  sheet.  The 
resulting  steady  state  RMS  wavefront  error  is  eomparable  to  the  eorreeted  wavefront 
results  aehieved  without  a  disturbance  but  at  the  cost  of  a  greater  settling  time,  Figure  61. 

The  filters  only  operate  well  if  the  frequency  of  the  disturbance  is  known  exactly, 
as  seen  in  all  the  notch  filter  experiments.  The  uncertainty  in  sampling  time  of  the 
experimental  testbed  prevented  the  notch  filters  from  performing  exactly  as  designed. 
Despite  the  filters  not  being  tuned  exactly  they  all  reduced  the  disturbance  regardless  if 
the  disturbance  signal  represented  actuator,  sensor,  or  struetural  disturbance.  The 
experimental  noteh  filtering  results  shown  in  Table  7  show  that  the  higher  bandwidth 
filters  tend  to  perform  better.  However,  a  wider  bandwidth  filter  may  have  some 
consequenees  in  a  non-laboratory  system,  where  adjacent  frequency  information  may  be 
desired. 
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Figure  60  Error  History  for  2.0  Hz  Sinusoidal  Disturbance  on  Wavefront  Focus,  with 
and  without  Notch  Filter,  Bandwidth  =  0.27i,  Initial  Bias  of  127 
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Figure  61  Error  History  for  2.0  Hz  Sinusoidal  Disturbance  on  Wavefront  Focus,  with 
and  without  Notch  Filter,  Bandwidth  =  O.Sti,  Initial  Bias  of  127 


Channel 

Frequency  (FIz) 

Bandwidth  of  NF 

o  without  Filtering 

o  with  Filtering 

10 

5 

0.571 

0.105 

0.074 

10 

2 

O.lTI 

0.224 

0.104 

10 

2 

0.571 

0.224 

0.075 

10 

2  and  0.8 

0.571 

0.465 

0.291 

10  and  16 

2  and  0.8 

0.571 

0.2185 

0.085 

10  and  11  (NF  on  10) 

2 

0.571 

0.312 

0.151 

10  and  11  (NF  on  10) 

2  and  0.8 

0.571 

0.329 

0.208 

10  and  11 

2  and  0.8 

0.571 

0.327 

0.134 

All 

2 

O.lTI 

1.778 

0.483 

Simulated  Vibration 

2 

0.271 

1.47 

0.708 

Simulated  Vibration 

2 

0.871 

1.47 

0.0405 

Table  7  Notch  Filtering  RMS  Wavefront  Error  Comparison 
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VII.  SUMMARY,  CONCLUSIONS,  AND  FUTURE  WORK 


A,  SUMMARY 

In  order  to  experimentally  validate  eurrent  teehnology  and  develop  new  wavefront 
eontrol  teehniques  for  future  spaee  telescopes,  an  adaptive  optics  testbed  was  designed 
and  built  at  the  Spacecraft  Research  and  Design  Center,  Naval  Postgraduate  School.  This 
testbed  is  unique  in  its  ability  to  experimentally  test  control  algorithms  for  both  jitter  and 
wavefront  control.  This  research  used  the  adaptive  optics  testbed  to  demonstrate 
wavefront  control  techniques  for  correcting  a  deformable  mirror  and  rejecting  known 
disturbances  introduced  into  the  mirror. 

Wavefront  estimation  and  representation  techniques  were  investigated  and 
evaluated  when  combined  with  an  iterative  feedback  control  algorithm.  The  control 
algorithms  were  evaluated  by  comparing  the  steady  state  RMS  wavefront  error  and  peak 
to  valley  wavefront  error  by  correcting  the  wavefront  from  a  biased  position.  The  control 
algorithms  were  then  evaluated  subject  to  a  5  Hz  sinusoidal  disturbance  on  the  actuators 
and  evaluated.  Notch  filtering  techniques  were  also  investigated,  with  known  sinusoidal 
disturbances  at  0.8  Hz  and  2  Hz  introduced  as  actuator  noise  and  as  a  simulated  structural 
disturbance.  The  notch  filter  was  combined  with  both  a  direct  and  indirect  iterative 
feedback  controller  and  the  RMS  wavefront  error  was  compared. 

B,  CONCLUSION 

This  research  has  shown  that  the  direct  modal  wavefront  estimation  techniques 
outperformed  the  indirect  wavefront  representation  techniques  when  combined  with  the 
iterative  feedback  control  laws.  The  indirect  wavefront  control  methods  could  be 
improved  by  using  a  singularity  robust  inverse  technique  when  inverting  the  influence 
matrix;  however,  the  direct  modal  wavefront  control  techniques  still  provided  superior 
performance. 

The  experimental  results  also  showed  that  a  combination  of  an  iterative  feedback 
and  gradient  feedback  control  law  using  variance  minimization  provided  the  smallest 
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RMS  wavefront  error  when  correcting  from  a  biased  position  and  when  correcting  the 
mirror  subjected  to  a  sinusoidal  noise  disturbance.  However,  the  combined  controller  is 
very  sensitive  to  the  gains,  requiring  both  the  iterative  gain  and  the  gradient  feedback 
gain  to  be  properly  tuned.  All  the  control  algorithms  showed  a  residual  disturbance  in  the 
wavefront  error  when  subjected  to  a  sinusoidal  disturbance. 

The  experimental  results  demonstrated  that  a  second  order  discrete  time  notch 
filter  can  be  used  in  the  adaptive  optics  control  algorithm  to  improve  the  steady  state 
RMS  wavefront  error  when  a  known  constant  frequency  disturbance  is  present.  Applying 
a  properly  tuned  notch  filter  in  series  with  either  a  direct  or  indirect  controller  decreased 
or  removed  the  dynamic  disturbance.  The  notch  filter  success  in  the  experimental  work 
was  dependent  on  knowledge  of  the  disturbance  frequency,  the  notch  fdter  bandwidth, 
and  the  knowledge  of  the  computer  sampling  time.  The  experimental  work  successfully 
demonstrated  that  low  frequency  filtering  of  actuator  noise  as  well  as  a  simulated 
structural  disturbance,  achieving  a  wavefront  error  comparable  to  results  achieved 
without  a  disturbance. 

C.  FUTURE  WORK 

The  current  testbed  and  future  research  can  be  improved  by  upgrading  the 
wavefront  sensors.  The  current  sensors  are  too  slow  to  observe  the  dynamics  of  the 
mirrors,  or  higher  frequency  disturbances.  Increasing  the  wavefront  sensor  speed  will  be 
an  immediate  improvement  to  current  and  future  control  algorithms. 

Although  the  control  algorithms  worked  well  on  the  current  testbed  the  current 
mirrors  do  not  have  the  dynamics  of  a  large  light  weight  mirror.  Future  research  should 
be  conducted  using  a  larger  light  weight  flexible  mirror  with  lower  structural  resonant 
frequencies.  This  will  require  robust  control  techniques  and  adaptive  fdtering  techniques 
to  control  the  system  when  the  disturbance  is  not  known  precisely,  and  to  account  for 
uncertainty  in  the  system.  Additionally  the  control  techniques  developed  in  this  thesis 
could  be  compared  with  multi-input  and  multi-output  control  techniques  that  incorporate 


80 


the  structural  dynamics  of  the  mirror.  In  addition  to  face  sheet  control  as  discussed  in 
this  paper,  particular  attention  should  be  given  to  segmented  mirror  control,  as  future 
space  telescopes  will  most  likely  be  segmented. 

The  SRDC  in  cooperation  with  the  Naval  Research  Labs  has  purchased  an  16  inch 
diameter  light  weight  six  segment  mirror.  This  mirror  will  allow  future  research  into 
control  techniques  to  properly  align  and  phase  the  mirror  segments.  In  addition  to  the 
new  mirror,  new  phase  diversity  sensors  will  also  be  included  in  the  experimental  testbed. 
Eventually  the  control  algorithms  discussed  in  this  research  may  be  combined  and 
adapted  to  be  used  with  the  segmented  mirror  control  algorithms  to  control  a  segmented 
mirror  with  deformable  face  sheets. 
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APPENDIX  A.  SOFTWARE  VERSIONS 


Software 

Version 

Funetion 

MATLAB 

R2007a 

Interfaee  Computer  with  MMDM, 
PDM 

MATLAB 

R2006a 

Interfaee  Computer  with  DSPACE 

DSPACE 

R5.1 

Interfaee  Computer  to  FSM 

Windows  XP 

Profession  Version  2002,  SP2 

Control  Computer  Operating  System 

Table  8  Software  Versions 
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APPENDIX  B.  SAMPLE  MATLAB  CODE 


A,  INITIALIZATION 

q, _ _ _ _ _  _  _  _ _ 

o 

%Matt  Allen  October  5,  2007 

%Iterative_inethod .  m 

%This  script  uses  Zernike  Polynomials  for  wavefront  estimation.  This 
%script  is  used  to  initialize  the  control  algorithm  by  computing  the 
%influence  matrix  (poke  matrix) . 

^  _  _______________  ___________ 

O 

clc; 

clear  all; 
close  all; 


%Define  Constants 


pw= .009630; 
fl=17. 361408; 
w2=wk; 


%pixel  width  in  mm 
%focal  length  in  mm 

%Weighting  Factors  for  Zernike  Polynomials  (Variance) 


%Load  and  process  Reference  Image  -  Calibrate 


TH=150;  %Threshhold  for  CCD 

I=imread (' Calibrate . bmp ') ;  %Read  in  calibrated  reference  image 

I ( 1 : 2 0 , 304 : 324 ) =0 ;  %Block  Out  Extra  Pixel 

[Grid  Ref  init]=Ref  Grid  init(I,TH);  %Compute  Reference  Grid 

Locations 

s_grid=Grid; 

%Ref erence_Grid=round ( s  grid);  %Round  Reference  Grid  to  nearest  pixel 
of f set_x=s_grid (64,1); 
of f set_y=s_grid (64,2)  ; 

s_grid ( : , 1 ) =s_grid ( : , 1 ) -of f set_x;  %Center  Reference  Grid  x=64 
s_grid ( : , 2 ) =s_grid ( : , 2 ) -of f set_y;  %Center  Reference  Grid  y=64 
Ref=s  grid; 


q, _ _ _ _ _ _ _ _ _ _ _  _ _ _  _ _ 

o 

%Build  Influence  (Poke)  Matrix 

q, _ 

o 

for  i=l:37  %Loop  based  on  Number  of  Actuators 

V=ones (1, 37) *sqrt ( . 5) *240;  %create  a  voltage  settings  vector 
V(i)=240; 

BAODMirror (hex2dec ( ' C800 ' ) , hex2dec ( ' CCOO ' ) ,  V)  ; 

pause (.1)  %Delay  Required  to  Command  Mirror 

BAOGrab(O);  %initialize 
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SH= (BAOGrab ( 1 ) ) ' ;  %grab  images  from  the  Shack-Hartmann 

WFS 

BAOGrab(2);  %de-initialize  mirror 

SH ( 1 : 2 0 , 304 : 324 ) =0 ;  %Blockout  extra  lenslet 

Data=Process  Image (SH, Ref  init)  ;  %Compute  Centroids,  pass  image, 
and  thresh  hold 

Data ( : , 1 ) =Data ( : , 1 ) -of f set_x; 

Data ( : , 2 ) =Data ( : , 2 ) -of f set_y; 

Zs ( : ,  i) =reshape (Data,  254, 1) ; 


end 


q, 

o 


%Estimate  Wavefront  of  Poke  Matrix  Using  Zernike  Polynomials  (21  terms) 
%return  Influence  matrix  of  Zernike  Coefficients 


B=Bmatrix ( Zs , Ref ) ;  %Builds  Modal  Influence  Matrix 

B.  WAVEFRONT  SENSOR  PROCESSING 

q, _ _ _ _ _ _ _ 

o 

%Matt  Allen  Aug  18,2007 

%Process_Image .m 

%This  function  is  used  in  conjunction  with  the  program  Ref_Grid_init .m 
%This  function  reads  in  the  reference  lenslet  location  'Ref_init' 
%computed 

%from  the  function  ' Ref_Grid_init . m ' and  sections  out  the  lenslet  array 
%areas  in  the  image  'SH'.  A  centroid  is  then  computed  and  returned  as 
%the  vector  'Grid'. 

^  _  __________  ________  ___  ___ 

O 

function  Grid  =  Process  Image (I, Ref  init) 

%Search  Sub-Image 
for  i=l:127 

P=double ( I (Ref  init (i, 2) -10 :Ref  init ( i , 2 ) +1 0 , Ref  init(i,l)- 
10:Ref  init (i, 1) +10) ) ;  %Section  out  Lenslet 
S=size(P);  S=S(1); 

M= ( sum ( sum ( P) ) ) ;  %Compute  Centroid 

X= ( linspace ( 1 , S , S ) ) ' ; 

Rx=P*X; 

Ry=X' *P; 

Rx=sum(Rx) *1/M; 

Ry=sum(Ry) *1/M; 

Grid (i, 1) =Rx-10+Ref  init(i,l); 

Grid (i, 2) =Ry-10+Ref  init (i, 2);  %X  and  Y  Locations  of 

Reference  Grid 
end 
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C.  CONTROL  ALGORITHM 


9-  _  ___  ___  __  __  _____  _____ 

o 

%Matt  Allen  31  October  2007 

%Control_Law.in 

0, 

o 

%Apply  Control  Law  Using  Iterative  Voltage  Updating  Vn+l=Vn-gBinv ( S ) , 
%where  S  is  the  slope  measurement.  This  algorithm  minimizes  the  slope 
%difference  as  directly  measured  from  the  Shack-Hartmann  Wavefront 
%Sensor.  The  algorithm  drives  measured  pixels  to  a  reference  pixel 
%associated  with  a  planar  wavefront.  This  script  requires  the 
%initialization  script  Iterative  method. m,  and  a  reference  image  of  a 
%planar  wavefront  (or  desired  wavefront)  named  Calibrate.bmp 

o  ___  ___  _ 

o 

clear  Error 
clear  Freq 

s=reshape (Ref, 254,1); 

BAOGrab ( 0 ) ; 

SH= (BAOGrab  (1) )  '  ; 

BAOGrab (2) ; 

SH(1:20,304:324)=0;  %Block  out  extra  lens let 

Data=Process  Image (SH, Ref  init) ; 

Data ( : ,  1 ) =Data ( : , 1 ) -of f set_x; 

Data ( : ,  2 ) =Data ( : , 2 ) -of f set_y; 

Data=reshape (Data, 254 , 1 ) ;  %SH  Location  Data 
SI=Data-s;  %Compute  the  Slope 

g= . 0 1 ; 

Binv=pinv(B) ; 

cold=Binv* (SI) ; 
tic; 

for  i=l:150 

BAOGrab (0);  %initialize 

SH= (BAOGrab ( 1 ) ) ';  %grab  image  from  Shack-Hartmann' 

BAOGrab (2) ; 

SH(1:20,304:324)=0;  %Block  out  extra  lens let 

Data=Process  Image (SH, Ref  init); 

Data ( : , 1 ) =Data ( : , 1 ) -of f set_x; 

Data ( : ,  2 ) =Data ( : , 2 ) -of f set_y; 

Data=reshape (Data, 254 , 1 ) ;  %SH  Location  Data 

SI=Data-s;  %Compute  the  Slope 

cnew=cold-g*Binv* (SI) ; 
for  p=l : 37 

if  cnew (p) <-l 


%initialize 

%grab  image  from  Shack-Hartmann' 
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cnew (p) =-l ; 
else  if  cnew(p)>l 
cnew (p) =1 ; 
end 

end 

end 

V=sqrt ( (cnew+1) * . 5) *240; 
cold=cnew; 

BAODMirror (hex2dec ( 'C800 ' ) , hex2dec ( 'CCOO '  )  ,  V'  ) 

pause ( . 03 ) 

end 
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